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I. INTRODUCTION 


A. BACKGROUND 

Surveying 1s the science and art of measurements which are necessary to 
determine the relative position of points above, on, or beneath the surface of the earth, 
or to establish the points in a specified position. Surveying operations are conducted 
not only on land, but also in the oceans and in space. The measurements of surveying 
consist of distances, horizontal and vertical, and directions. In order to provide a 
framework of survey points whose horizontal and vertical positions are accurately 
known, basic horizontal and vertical control surveys must be performed. A primary 
use of control surveys 1s for construction of control for a map or chart. The 
fundamental network of points whose horizontal positions have been accurately 
determined ts called horizontal control (Schmidt, 1978, p. 122]. 

Horizontal control generally is established either by traverse, triangulation, or 
trilateration. Which one 1s to be used depends on the accuracy required and the factor 
of economy in the selection of survey method. Obviously, there are many degrees of | 
precision possible in any measurement because no surveying measurement 1s exact. 
Each of these methods may be the best one to use for a given purpose. Ordinarily, it 1s 
a waste of time and money to attain unnecessarily high accuracy. On the other hand, 
if the measurements are not sufficiently precise, faulty survey results are produced. 
Therefore, the best surveyor is not the one who makes the most precise measurements, 
but the one who 1s able to choose and apply the appropriate measurement with 
precision requisite to the purpose. 

Before 1950 the main framework of a first-order geodetic survey almost always 
consisted of triangulation, which could be replaced by traverse in cases where the 
topography made triangulation impracticable. Today, due to the development of 
electronic distance measuring (EDM) equipment, the first-order control points can be 
established by means of high accuracy traverse [Allan, 1968, p. 370]. Therefore, the 
horizontal control is frequently provided by traverse, especially for surveys in area of 
limited extent, mostly flat and jungle covered. Traverse in such cases is much more 
economic, convenient, and rapid than other methods and the results are equally 


acCCurale. 


In order to achieve high precision of horizontal control points when distributed 
over a large area, first or second-order geodetic surveys are required. These types of 
survey treat the shape of the earth as ellipsoidal and require using the most accurate 
distance and angle measurements. Computation of such a survey is relatively 
complicated, based on long geodetic formulas for computing (with necessary precision) 
the exact horizontal and vertical position of widely distributed points on the earth’s 
surface. 

Disregarding ellipsoid shape, a third-order survey is used over earth areas of 
limited extent. In this type of a survey, the earth can be considered as flat and all 
angles are considered to be plane angles. Surveys of this type are used in the 


densification of geodetic control. 


B. OBJECTIVES 

As mentioned above, the traverse method has been used worldwide mostly for 
densification of control stations. However, there are many methods of traverse 
computations. The main objectives of this thesis are to (1) compute a 
closed-connecting traverse and adjust station positions by the Approximation Method . 
with the compass (Bowditch) rule and Least Squares Method (adjustment by 
observation equations only), and (2) compare the results of the two methods. 

All computations have been accomplished in the Universal Transverse Mercator 


(UTM) grid coordinates rather than geodetic coordinates. 


I. TRAVERSES 


A. GENERAL 

The word traverse generally means to pass across. But tm surveying, this word 
means the measurement in a specified sequence of the lengths and directions of lines 
between points on the earth whose position may be know or unknown. Traverse ts the 
most Widely employed method for densification of local horizontal control. Linear 
measurements are made either by direct observation using a tape or EDM cquipment, 
or by indirect observation using tacheometric methods. The angular measurements are 
made with theodolite or transit. In this thesis, the only traverse operations considered 
are for angles measured bv theodolite and distances measured by precise 1: DM 
equipment or tapes. 

Two kinds of traverse exist in surveying, an open and a closed traverse. 

1. Open Traverse 

An open traverse normally ortginates at a pomt of known position and does 

not return to the starting point nor does it termmate on another point of known 
position (igure 2.1). Open traverses should generally not be used because they can 


not be checked for errors. 


Known station 


unknown station 


measured distance 


observed angle 





Figure 2.1 -Open: liravense: 
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2. Closed Traverse 

Closed traverses can further be sub-divided as, a closed-loop and a 
closed-connecting traverse. 

A closed-loop traverse is one that originates and terminates on a single point 
of known position, thus forming a closed polygon (Migure 2.2). Vhis type of traverse 
provides an internal check on angles but no check on systematic errors in distance. 
Also, if the starting azimuth (between stations 0 and [ in Figure 2.2 ) has an error, it 
causes crror in orientation of the entire traverse. A closed-loop traverse generally 


should not be used. 





A 6 = known station 
© = unknown station 
in; = measured distance 


fi = observed angle 


Frgure 2.2) Closed-loop Traverse. 


A closed-connecting traverse ts one that begins and ends on two diflerent 
points whose horizontal positions have been previously determined by a survey of 
higher or equal accuracy (Figure 2.3). ‘This type of traverse is preferable to all others, 
since computational checks are possible to detect systematic errors in both distance 


and direction. 


B. ANGULAR AND DIRECTIONAL MEASUREMENTS 

The position of traverse points is determined by the direction and distance from 
the starting point. To obtain the direction by means of azimuth, the horizontal angle, 
or plane angle, must be measured in the field. Also, the deterinination of vertical 
angles, or zenith distances, may be required to reduce slope distances to horizontal 


distances. 


= known station 


unknown station 






2 ©) 


= measured distance 


Ai = observed angle 


ligure 2.3 Closed-conmeetimomlinan erse: 


Commonly, both horizontal and vertical angle measurements are accomplished 
with a transit or theodolite. The theodolite is employed especially for surveys of high 
precision. Two types of theodolite are a repeating theodolite and direction theodolite.  - 
A repeating theodolite reads directly to 20” or I’ and may be estimated to one-tenth of 
the corresponding direct reading. A direction theodolite usually reads directly to 1” : 
and may be estimated to 0.1” [Davis et al., 1981, p. 215]. In general, a direction 
theodolite is more precise than a repeating theodolite and with it, plane angles are 
computed by subtracting one direction from another. 

In all types of traverses, the horizontal angles can be measured by one or more 
of the following listed angle measurement methods. 

1. Interior Angle 

Interior angle is the angle measured within a closed figure at the intersection 
of two lines (I‘igure 2.4). 
2. Deflection Angle 
Deflection angle is the angle measured from the extension of the preceding line 
to the succeeding line (figure 2.5). Such angles must be identified as right or left to 
express whether the angle is turned to the right or to the left from the preceeding line. 


3. Angle To The Right 


Angle to the right is the clockwise angle measured from the preceding line to 


the succeeding line (Iigure 2.6) 


2 = observed interi 
acilie rior 





Figure 2.4 Measuring of Interior Angles. 





Figure 2.5 Measuring of Deflection Angles. 


C. LINEAR MEASUREMENT 
Direct lincar measurements may be obtained in traversing by pacing, odomcter 
reading, tachcometry (stadia), subtense bar, taping, and EDM. Of these methods, 
taping and EDM are most commonly uscd by surveyors. !lowever, [IDM equipment 
has a clear superiority over traditional taping for lines in excess of about 250 meters. 
[Distances measured using EDM cquipment are subject to personal, instrumental, 


and natural crrors. Personal crrors include misreading, improper centering of the 





igure 2.6 Measuring of Angles to the Right. 


mstrument over the stations, failing to exactly center the null meter, and incorrectly 
measuring meteorological factors and instrument herghts. Instrumental crrors, 
expressed in terms of the accuracy of the mstrument specified by the manufacturer, : 
contain two parts. For example, if the accuracy of an instrument ts designated as + 

(10 ppm + 5 mm), the constant error part is + 5 mm, which ts independent of the 
distance, and the value of the proportional part is [0 ppm (parts per million) which ts a 
function of the distance measured. Constant error is most significant for short 
distances. lor very long distances the constant error becomes negligible, but the 
proportional part 1s important. Natural errors such as refraction are cause by 


changmg of atmospheric conditions along the measured line between the end stations. 


D. ACCURACY 

In survey adjustment, a deviation from the ‘true’ value ts considered as an 
observational error and the standard error designates the measure of accuracy of the 
observation. The meaning of an accuracy ts then the degree of conformity or closeness 
of a measurement to true value. 

The quality of traverse operations ts dependent upon the accuracy of angular and 
linear measurements; thus, in checking the accuracy of traverse two quantities are 
considered, the angular misclosure and the linear misclosure. Although the positional 
closure (relative accuracy) 1s an indication of the overall quality of the traverse and ts 
used for traverse classification, it does not yield information on the precision of point 


location determined in a traverse [Davis et al., p. 332]. 


Id 


The inherent weakness in a traverse is that the deviation of each measured line is 
determined by a single series of angular observations, further, any error in any angle 
will affect not only the adjoining line but all subsequent lines to a greater or lesser 
extent according to their lengths [Allan, 1968, p. 371]. 

Angular misclosure 1s expressed by standard error of the measured angle times 
the square root of the number of measured angles. 

Linear misclosure is commonly expressed as a ratio of total misclosure to total 
length of traverse. 

Finally, some of the most significant features of traverse classification by the U.S. 


Federal Geodetic Control Committee (1984) are shown in Table I. 


FE. ADJUSTMENTS 
Adjustment ofa traverse is carried out to ensure consistency within the known 
positions of the originating and terminating stations and to remove inconsistencies in 
observed angles and distances to compensate for random errors. For a more precise 
extended traverse, adjustments made on the basis of least squares are preferred. Buta - 
traverse of limited extent can be adjusted by simple approximation methods. 
1. Approximation Methods 
There are four methods for traverse adjustment by approximation. 
a. Arbitary Method 
This method does not conform to a fixed rule. Rather, the linear error of 
closure is distributed arbitarily according to arbitary preference of the surveyor. 
b. Transit Rule 
Transit rule 1s better for adjustment of the traverse where the angles are 
measured with greater accuracy than distances, and 1s valid only when the traverse lines 
are parallel with the grid system used for the traverse computations. Corrections are 
made by the following rules: the correction in latitude for any station 1s equal to the 
multiple of latitude in that section and total closure in latitude divided by the sum of 
all latitudes in traverse, and the correction in departure is equal to the multiple of 
departure in that section and total closure in departure divided by the sum of all 
departiimes netiertraverse (Mavis ct al., 1981, p. 323]. 
c. Compass or Bowditch Rule 
This method 1s suitable for adjustment of the traverse where the angles and 


distances are measured with equal precision and uses the following rules: the correction 


Abie 
TRAVERSE CLASSIFICATION 


Order First Second Second Third 
Class I ee i 


StaAClOn y= bactad solr 
less than (km) 10 & Z 


Minimum number of 
NETWOrk ConcroL 
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Theodolite least 
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Azimuth closure at 
azimurch check worn 
(second of Sins) 
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after azimuth One On OF Os oie 
adjustment* le JOOO00O0O 1.50000 Te 20000 aatoOoCcor!: 50050 
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MEQUAREE ANCL se keer Coe ct ed eae 

ae closure (e.g., 1:50,000) is relative error of 
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in latitude for any station is equal to the multiple of the length in that section and total 
closure in latitude divided by the total length in the traverse, and the correction in 
departure is equal to the multiple of the length in that section and total closure in 
departure divided by the the total length in the traverse [Schmidt, 1978, p. 150]. 
d. Crandall Method 
Crandall method is a rather complicated procedure which is more rigorous 
than either the compass or transit rules but suitable for adjusting traverses where the 


linear measurements contain larger random errors than the angular measurement. In 
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this method, the angular error of closure is first distributed in equal portions to all of 
the measured angles, then linear measurements are adjusted by using a weighted least 
Squares procedure [Brinker, 1977, p. 228]. 
2. Adjustment by Least Squares Method 

The method of least squares adjustment is based upon the theory of 
probability; it simultaneously adjusts the angular and linear measurements to make the 
sum of the square of the residuals (error) a minimum [Brinker, 1977, p. 228]. This 
method can be used for any type of traverse. Because of the availibility of fast 
computing devices at the present time, the least squares method is being widely used. 
Further, the least squares solution has the advantage that it determines, quite 
objectively, a unique solution for a given adjustment problem [Clark, 1973, p. 121]. 

In general, adjustment 1s needed whenever there are redundant observations 
(more observations than are necessary to solve the required unknowns). As an 
example, to determine the angles of a plane triangle, only two observed angles are 
required because the third angle can be obtained by subtraction from 180°. When 
three angles are observed, the sum of them will not be equal 180° due to error in 
measurements. Therefore, these three angles should be adjusted to fit the functional 
model. 

The redundancy may be interpreted to mean that among n observations there 
exist r conditions or functions (n > r) that must satisfy the model. 

Let n be a number of observations and no the minimum number of 
observations to find the uniquely solution in the model, then redundancy or degree of 


freedom in the statistic, r, 1S 
ce (23) 


Consequently, there are r redundant observations, which can also give a 
solution. To detect the error in each observation, the best estimated or the most 
probable value must be defined because the true value is not known exactly. 
Statistically, the best estimated value of a group of repeated observation 1s the average 
(arithmetic mean). 

Once the difference between observed value (X,) and estimated value (X,) is 
determined, the adjusted value (X,,) 1s obtained through a least squares solution, then 


the residual (v) can be expressed as 


ee 


v= X,-X (272) 
and 
— = se steal (233) 


where dx is the correction to estimated value to obtain the adjusted value. 

The least squares adjustment method is based upon the criterion of the sum of 
the squares of the observational residuals must be minimum. 

When observations are considered as uncorrelated and of equal precision (with 


identity weight matrix), the least squares condition can be expressed as 

MD = ak ea = Wc = minimum (2.4) 
Or in matrix form 

® = V!v = minimum (2.5) 


where V is the vector of residuals. 
In uncorrelated observations with unequal precision, such as distances and 
angles [Mikhail, 1981, p. 68], the Equations 2.4 and 2.5 become 


MD = WV) + Wov97 +. FW hve = Vw. = minimum (2.6) 
Or in matrix form 


MD = viwv = minimum (235 


where w; is the it element of the diagonal weight matrix W and v, is the residual 


1 
associated with the corresponding i") observation. 

Generally, the relative weights are inversely proportional to variance, thus the 
weight matrix is the inverse of cofactor matrix, Q (when it is square and nonsingular) 


and defined as 


Ww = Q! (2.8) 
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where the elements of cofactor matrix Q are 


Qi = 9;7/69" (2.9) 
and 
ai = G;,/50° (2.10) 
y th 


where oF the variance of the 1°" observation, oi; is the covariance between the ith and 
jth observations, and os is variance of unit weight [Mikhail, 1981, p. 67]. 

For the case of uncorrelated weight observations, the cofactor matrice will be 
diagonal with all off-diagonal elements being equal to zero, thus the diagonal elements 


of weight matrix in this case are 


ee eg = Sy7/s.. Parl) 

Generally, there are two types of equation 1n least squares adjustment: 
condition and observation equations. Condition equations include one or more 
observations but observation equations include parameters and only one observation. 

The condition as well as the observation equations involved in an adjustment 

problem can be linear or nonlinear. However, least squares treatments are generally 
performed with linear functions, since it 1s rather difficult and often impractical to solve 
nonlinear models [Mikhail, 1978, p. 108]. Consequently, whenever the equations in the 
model are originally nonlinear, they have to be linearized. A method of series 
expansion, especially Taylor’s series, is often used to obtain linear equations. Only the 
zero and first-order terms are used and all other higher-order terms are neglected. 


Thus, the linearized form for the general case of m functions of n variables becomes 


— FO 
= a eK (2.12) 


where F® is the zero-order terms, when x = x0 at ee a is a Jacobian matrix of 
coefficients of first order of n variable (Appendix A). 
The choice between the observation equations (indirect observation) and 


condition equations (observation only) techniques will depend mainly on the 


19 


mathematical model of the problem to be solved. However, the final answers are 
always the same. - 
a. Adjustment by Observation Equations 
The method of adjustment by observation equations is performed with both 

observations and parameters. Therefore, the number of equations is equal to the 
number of observations. Using the example at the beginning of this section, if three 
measured angles and their residuals in a plane triangle are assummed to be @, B, y, cae 
V9, and oe respectively, and the adjusted value of those angles is ey and Xn3 
then the three condition equations in the normal form (Zero at the right-hand side) are 

Ce ie | 

(Paava) Xe oe 0 

Cea) oa = 0 
by using Equation 2.3 for adjusted values 

V,- Ope X.j- @ 

ee aca 

oe 























letting 
V -1 O O ax x - 
v= l|vdi!, B= | 0-1 O], x = /axd|, £ = |xel - B 
vs Oona dx ez _ 
3 3 Xo37 7 


then, expressed in matrix notation as 

V3, + B33X3, = F3, and the general form of the adjustment of observation 
equations become 
Ep 


Vii +BayX (2.13) 


Wl nl 


where V is an n by | vector of residuals; B is an n by u matrix of coefficients; X is an u 
by 1 unknown vector in which u is the number of parameters; and fis an n by 1 vector 
and equals P - X., in which P is evaluated using the estimated value X,. 
b. Adjustment by Condition Equations 
The method of adjustment of condition equations only has no parameters 
included in the condition equations. Thus, the number of condition equations is equal 
to the number of redundancies. From the example above, n= 3, Ng Zea Te 
which is the number of condition equations being set. In this case, because the sum of 
interior angles must equal 180°, the single condition equation can be expressed as 
ET Dismay PN a = FS0e 
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Vi tVvotv, = 180°-a-f-y 


5 Vv 

v2 
ec’ 
Then, the general form of this technique 1s 


oe aly V= i F = |180-a-B-y| 








“i (2.14) 


When the conditions are originally linear, the vector F is usually written in terms of the 


given observations as 
ee aus tno! (2.15) 
where A is the coefficient matrix V., P is a constant term (see Section II.E.a), X, is 


observed values, r is redundancy, and n is a number of observation [Mikhail, 1976, p. 


Ly: 
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II. DATA ACQUISITION AND REDUCTION 


A. DATA ACQUISITION 

Taverse data used in this thesis were obtained from field work accomplished from 
25 September thru 9 October 1972 by CAPT Glen R- Schaefer, NOAA Conps) ancien: 
Jim D. Shea, National Ocean Service (NOS), utilizing traverse methods in Pinellas 
County, Florida. Only the first 15 of 40 occupied stations and three intersection points 
will be used for analysis (Figure 3.1). The two pairs of known stations for this 
closed-connecting traverse are shown in Table II. 

The known stations were observed by the US Coast and Geodetic Survey (now 
NOS) and adjusted by the National Geodetic Survey (NGS). Station Turtle 2 1s of 
first-order and the other three stations are of third-order. 

The horizontal angles measured by the method of angles to the right (Table ITI) 
Were turned with a Wild T-2 theodolite, according to the specifications of third-order 
class I traverse, by starting at station Tomlinson and using Egmont Key Lt. House for , 
a backsight. The traverse was closed on Turtle 2 with a check azimuth to Madeira 
Beach Tank. 

The slope distances were measured in the field with a Model 76 Geodimeter in 
feet and corrected for temperature and pressure. Distance measurement by the Model 
76 Geodimeter are reported to have an accuracy in the temperature range of -20°C to 
+ 50°C of £(1 ppm + | cm) with a resolution of 1 mm. [Schmidt, 1978, p. 116]. all 
observed distances were converted to meters and reduced to horizontal by the 
procedure given later in this chapter. Finally, geodetic distances are reduced to grid 


distances by applying the scale factor correction (Table IV). 


B. DATA REDUCTION 
1. Computation of Ellipsoidal Distances 
For the requirement of high precision in the first-order traverse, the measured 
distances obtained by EDM equipment are first corrected for atmospheric conditions 


and then reduced to the ellipsoid (Figure 3.2) by the equations 


S = 2R, Sin )(Sp/2Ry) (3.1) 


Ze 
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Figure 3.1 Traverse Layout. 
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TABLE II 
DATA OF KNOWN STATIONS 


Station Station 
No. Name 


O Egmont Key Lt. House 
1 Tomlinson 

> OE eer 

6 Madeira Beach Tank 


Hl 
1 





Sy = ([m?-(hy-hy)] / (1+ hy/Rg (1+ h/Rg) }1/? (3.2) 
I/Rg = {(Cos*a)/M} + {(Sin2a)/N} (3.3) 
M = a(I-e”) / (1-e7Sin2@)>/2 (3.4) 
N = a/(1-e7Sin2g)!/2 (3) 


where m is the measured distance, S is the ellipsoidal distance, Sp is chord distance, hy 
and hy are the heights above the ellipsoid at each end of the line, Rg is the radius of 
curvature of the chord distance m at an azimuth @, @ is geodetic azimuth clockwise 
from north, M is radius of curvature in the plane of the meridian, N is radius of 
curvature in the prime vertical, a is a semimajor axis of the reference ellipsoid, e is its 
eccentricity, and @ is the mean latitude of that line. [Torge, 1980, pp. 48,49,50,51,180] 
2. Computation of Geodetic Distances 

For a lower-order traverse, the measured distance can be reduced to mean sea 

level (MSL) or geoid only. Because the error of only 1 ppm in length results from an 


error of 6 m in separation between MSL and the spheroid [Bomford, 1980, pp. 42, 342, 
345]. 


TABLES 
OBSERVED HORIZONTAL ANGLES 


A) For Traverse 


At Stn. Stn. ge 0 
O. Name D S (Ce: ) 
QO Egmont Key Lt. House 
1 Tom! Lnsen 105 21 25 5" 
2 TN-O1 243 39 18 5" 
3 TN-O02 L6oCuewerOn 15 5" 
4 TN-0O3 Somes) 56 5" 
5 Ruscue 291 28 29 ae 
6 RE-Ol MoOwnes, 42 5" 
7 RE-0O2 Zeo9F 55 50 5" 
8 RE-03 92 43 29 5" 
9 RE-04 yee. 51 5" 
10 RE-05 182 31 19 5" 
A. RE-06 196 54 42 5" 
LZ RE-08 168 51 46 5! 
3 RE-0O9 161) 0G 42 Be 
14 RE=-10 Zoom oom 14 ae 
iS Due eee 78 #37 24 5" 
16 Madeira Beach Tank 
Eye Onrelmcernseccion Points (IP) 
IP #1 * 
At Traverse endl ee CO 
Back Stn. Sel. D S (+) 
RE-O1 RE=-02 196 23 46 By! 
RE-04 RE=-05 hod S40) 22 Bh 
RE-0O5 RE-06 T3627 oe 5 5" 
RE-06 RE-08 Sa 45eu03l1 5! 
RE-08 RE-09 Zoe eee 2 2 5" 
LEA 
RE-03 RE-04 nS S57 Ol 5" 
RE-04 RE-05 97-12 06 5" 
RE-05 RE-06 Ii SI MAE 5" 
IP oO kKk* 
RE-08 RE-09 N94 Ol) 22 5" 
RE-O9 RE-10 237 ee). oe a 


See ecersourg SCH CO Tank 
eee cer obuUrd —bGH Sty Jeans CH Tower 
xx* Bay Pines Veterans Administration Hosp. 
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m = measured distance m 

So = chord distance h, 

S = are distance on h, 
ellipsoid 

h, = elevation above 


ellipsoid at 
Station 1 


h, = elevation above 
ellipsoid at station 2 


Ro = radius earth curvature 
along measured line 0 


Figure 3.2 Lllipsoidal Distance. 


lfowever, the process of reduction requires three steps: (1) correct slope 
distances to horizontal distances, (2) reduce horizontal distance to geodetic distances, 
and (3) change the geodetic distances to grid distances. 


The slope distance data used included correction for atmospheric conditions. 


In a plane survey, such as a traverse, there are two considerations for the 
reduction of slope distances to horizontal distances, a short slope distance and a long 
Slope distance. 

a. Slope Reduction for Short Distances 

Short slope distances (S 2 mi or 3.3 km) measured by using EDM 


instruments separate from the theodolite, can be reduced to horizontal with a simple 


trigonometry process as 
D = imcosQ@ (3.6) 


where D is the horizontal distance, m ts the slope distance, and 0 is the vertical angle 
(Figure 3.3). 
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‘slope correction 
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Figure 3.3. Slope Reduction for Typical Triangle. 


The horizontal distance also may be determined by using the difference in 


elevation between the two ends of the line. The horizontal distance ts 
D = (m*- d2)!/2 (3.7) 
in Which d is the difference in elevation between the two end points. The heights of the 


EDM instrument and reflector above the survey mark must be observed, and d 


becomes 


oe 


d = (E,+1,)-(E5+0,) (3.8) 


where E, and E, are the elevation at each end of the line respectively, I, is the height 
of the instrument, and O, is the height of the reflector. Then, expanding the right side 


of equation 3.7 with the binomial theorem yields 
D = m-(d2/28 + d4/8S? + ...) (3.9) 


The quantity in the parenthesis is designated slope correction. For 
moderate slope the first term is usually adequate. When the slope distances and 
vertical angles are obtained by separated EDM equipment from theodolite, the 


correction of the vertical angle must be determined. The corrected vertical angle 0+ 1S 
0, = O54 Ad (3.10) 
where @p is an observed vertical angle ty theodolite, and A6 is 
A9 = H Cos @)/S Sin 1” (37108 
and 
PAs = (Aa aol ae ont (3:12) 


where H, is the height of the reflector, H, is the height of the target, H,; is the height of 
the EDM, and H, is the height of the theodolite [Davis et al., 1981, pp. 103-104]. 
Equations 3.11 and 3.12 are not needed when the slope distances and vertical angles 
are obtained simultaneously by using an EDM transmitter built into a theodolite. 

b. Slope Reduction for Long Distances 

Slope reduction for long distances (> 2 mi or 3.3 km) involves using 

vertical angles affected by curvature and refraction. By assuming a mean radius for the 
earth of 3959 mi or 6371 km, then the curvature correction (C), expressed as an angle 


in seconds, 1s 4.935” per 1000 ft cr 16.19” per km and the horizontal distance, D, is 


D = m Sin(90°-0-C) / Sin(90° + C) (Gag) 
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for 
= (yp) 2 (3.14) 


where y and f are the vertical angles at each end of the measured line [Davis et al., pp. 
106-107]. 

When a single vertical angle (y) is observed, @ is the corrected vertical angle 
for combined results of curvature and refraction (C&R), then, 8 is y + (C&R). The 
C&R correction is 4.244” per 1000 ft or 13.925” per 1000 m. The correction of C&R 
will be positive when the vertical angle is an elevation angle and negative in the case of 
a depression angle [Davis et al, 1981, p, 108]. 

3. Reduction of Horizontal Distances to Geodetic Distances 
The horizontal distance at same elevation above the geoid, must be reduced to 


a geodetic distance. This can be done by the equation 
D> (8D) (Rte) (3.15) 


where D’ ts the geodetic distance, R is the mean radius of the earth’s surface at that 
section, D is the horizontal distance at elevation E above the geoid [Davis et al., 1981, 
ore Oar 


C. GRID DISTANCES 

The traverse computation, based on the UTM grid coordinate system, requires 
the reduction of geodetic distances to the plane of the projection by applying the 
projection scale factor and grid scale constant. Scale factor can be obtained from a 
graph or from a rigorous formula [Department of the Army Technical Manual, 1958, 
Bp Oc lae 1.€:, 


lee < VIIT) q? + 0.00003 q? | (oa) 
where k is the scale factor at scale working on the projection, Kg is the central scale 
factor which is an arbitary reduction applied to all geodetic lengths to reduce the 


maximum scale distortion of the projection (for UTM, kg = 0.9996), and values for q 
and (XVIII) are obtained by the formulas which shown in Table V. 


IB 


TABLE V 
SPECIPICA TION OF PAs | iS 


ice er Cos“@ : . 10% 
Z v2 
a 
(1 - e* Sin*g)1/2 
(ecentricity)? = (ar - b“) vy at 
ef / (1 - e*) 


© COCOCHS, 


grid easting = E' + 500,000 when point is 
east of Central meridian, 30073000. 
when point is west of central meridian 


radius of curvature in the prime vertical 
semi-major axis of the ellipsoid 
semi-minor axis of the ellipsoid 

latitude 
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IV. TRAVERSE COMPUTATION AND ADJUSTMENT 


Linear measurements and angles must be checked by computation to determine 
the position of traverse stations and whether the traverse meets required precision. 
Traverse station coordinates are usually expressed in terms of geographical coordinates 
Or rectangular coordinates such as those based on an UTM projection. Traverse 
computations are usually done in rectangular coordinates because of the ease of 
computation. In this thesis, only closed-connecting traverse computations in UTM 
coordinates were used. When specifications were satisfied, the traverse was adjusted 


for perfect geometric consistency among angles and lengths. 


A. DATA PROCESSING 
1. Set up of data base 


Two files on IBM 370/3033AP main frame computer at NPS were established. - 


2. Modification of an Existing Program 
The Fortran program TRAVAD,, originally written by LCDR Saman 
Aumchantr, RTN (1984) in Watfiv language for computing and adjusting the traverse 
station coordinates, was modified to be able to handle the distances reduction 
processes. 
3. Writing a New Program 
A Fortran program INDTRA was written for computing and simultaneously 
adjusting traverse station and intersection point coordinates by least squares 


observation equations method. 


B. COMPUTATION OF STARTING AND CLOSING AZIMUTHS 

The directions of lines by means of azimuth are used for traverse computation, 
because sines and cosines of azimuth angles automatically provide correct algebraic 
signs for latitudes and departures. 

The terms latitude and departure are widely used in rectangular coordinate 
calculations of surveying. The latitude of a line is its projection on the reference 
meridian, which differs from geographic latitude. The departure of a line 1s its 
projection on the east-west line perpendicular to the reference meridian. In traverse 
calculations, north latitudes and east departures are considered plus; south latitudes 


and west departure, minus. 


Bil 


Latitudes are also sometimes termed ‘northings’ and ‘Y differences’ (AY); 
departures are simuarly called ‘eastings’ and *X differences’ (AX). 

Because the closure angle of traverse can be checked by azimuth of each 
consecutive line, starting and closing azimuths have to be determined in the first step 
of computation and azimuths can be computed from a pair of Known coordinate 
station positions at the two ends of the traverse (Figure 3.1) 

The azimuth of the line from A to B, QAP: measured clockwise from north, is 


determined by the equation 


Gap = Tan (AX / AY) (4.1) 
for 

eon (4.2) 
and 

AY = Yp- Ya (4.3) 


where Xq and Xp are the grid easting coordinates, and Yq and Yp are the grid 
northing coordinates of stations A and B, respectively (Figure 4.1). The quadrant of 
the azimuth of line AB, Gap, is dependent upon the sign of AX and AY (Table VI). 
The back azimuth Gp q (the azimuth from B to A) is obtained by adding 180° to the 
forward azimuth Gy p. 

The length of the line AB (denoted as Sap or S) can be determined by the 


Pythagorean theorem or by one of the trigonometric relationships 


nn 
I 


AX / Sin a (4.4) 


Or 


nN 
I 


AY / Cos @ (4.5) 
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Figure 4.1 Aztmuth Computation. 


TABLE VI 
CUADIOVNT OF AZ tl 


Quadrant of azimuth measured Sign Of 
Clockwise from north AX 


2 
8 
7 
6 


O 
180 
270 
300 





Substituting data from Table I into Equations 4.2 and 4.3, the AX and AY 
between Egmont Key Lt. House and Tomlinson are computed as AX = 3185.585 m 
and AY = 9470.704 m. The azimuth of the line from Egmont Key Lt. [louse to 
Tomlinson is 18° 35’ 27.6” (Equation 4.1). Similarly, the azimuth from Turtle 2 to 
Madeira Beach Tank is computed to be 273° 57° 12.0". These starting and closing 
azunuths will be used for computing the coordinates of cach traverse station and for 


checking the angular error. 


Sp 


C. COMPUTATION OF TRAVERSE STATION COORDINATES 

Computation of traverse station coordinates is the reverse process of finding 
azimuth and distance from coordinates. Therefore, the rectangular coordinates for 
each closed-connecting traverse station can not be computed unless forward azimuth 
and distance from the previous station are known. 


The azimuth is reckoned clockwise from north and obtained by 


where ik is the forward azimuth from station j to station k, ai, is the forward azimuth 
from station 1 to station J, and B; is the horizontal angle at station j for j values of 1 to 
n, where the previous station i = j - | and the next station k = j + 1. The number } 
will increase from | (which designates the starting known station of the traverse) to 
number n (which was the last occupied and known station of the traverse). 

Departures and latitudes are then computed by using Equations 4.4 and 4.5 


which are rewritten as 


and 


where Sik is the distance between stations } and k. The coordinates of all other 
traverse stations can be determined by adding successive departures (AX5,) and 
latitudes (AY) to the previous station’s X and Y coordinates, respectivelv. 

Using the data in Tables II, III, and IV, the azimuth and coordinate 
computations at the first station are shown here by using Equations 4.6, 4.7, and 4.8. 


The result for other stations can be seen in Table VII. 


Computing the azimuth from station 1 to 2 
AZ = 18°35°27.6 + 180° + 105 21 25 = 8303 seer oe 


Computing the coordinates at station 2 
X-easting = 329400.420+ [701.318 Sin(303°56’52.6”)] 
= 328818.645 m 
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Y-northing = 3063485.483 + [701.318 Cos(303°56’52.6”)] 
~ = 3063877.126 m 


When the coordinates of all stations have been computed, they are still 
unadjusted coordinates and can then be adjusted by one of the two methods of Section 


I{.E (approximation and/or least squares method). 


TABLE VII 
UNADJUSTED TRAVERSE STATION POSITIONS 


stn. Angles Forward ist. Unadjusted 
Azimuths Coordinates 
baetie <S DaeeM Ss X(m) Y(m) 


Mey Sly MAS) 
ces 2i.°25 *329400. 420 *3063485. 483 
8105 “56%53 e318 
243-39 18 S206 1645 3063877. 12Z6 
yoo 1 4.090 
HTS cl eal 328923.709 3064664. 236 
355 46 26 . 068 
Sooo SO Bige4 OO. 30660017. 613 
ZSSe4222 Ze 
Zoe Ao 29 S20 toMmOZ9 soUGsS95. /60 
347 10 51 SOF 
160 43 42 Solero 3007570. O99 
OZ) 34°35 2c 
209. 55> '50 SZueoluasoe BCO9106. 259 
5/ er 23 . 481 
227 ao 29 OZ 7iooS oe 5009374. 789 
33Q@33 52 . 448 
7S 2265 327418.001 3069801. 054 
S265 50 205 oS 
ez 31 19 ey Cops ZS O070353..210 
6 Peo az : 
196 54 42 S20 /+sel6 9530/0985. 723 
348 22 44 oo 1 
168 51 46 SZC lS moomoo 71391. 785 
So 7512 30 - 415 
161 6 42 SZe061. 476 93072814. 113 
Sie o2 be Zo 
236 58 14 tole oto 6 50/4219. 797 
oo 26 346 
Lomo 2 & SVoc ce wloO wm 30/61/8924 
273) 50, 50 


0 
a 
Z 
3 
5 
6 
i 
8 
9g 


RPM Be F RP 
mm bp W NYO # OC 


+ 
OV 


Coordinates for station 1 were known and held fixed. 
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D. ADJUSTMENT BY APPROXIMATION METHOD 

In this thesis, the method of Compass or Bowditch rule was used to adjusted the 
data in Tables III and IV. Thus, the first step 1s to determine the angular error of 
closure and adjust the angles to obtain the proper closing azimuth (closed azimuth at 
fixed stations). 

1. Angular Errors of Closure 

In the closed-connecting traverse, when there are n stations of observed 
horizontal angles, n-! lines will be measured. An angular error in traverse can be 
checked and obtained at the last station by comparing the computed azimuth and 
closing azimuth at the known station. 

The closing azimuth computed (from the Known station coordinates at the 
traverse end) at the station [5 1s 273° 57’ 12.0”. But because of error in measurement, 
the azimuth computed through the traverse at this station is 273° 56’ 49.6”, which is a 
difficiency of 22.4”. This amount of angular error in 15 observed stations meets the 
limit for allowable error for a third-order class I traverse (allowable error from Table | 
is 38.7”), 

The average correction (Table VIII) is distributed uniformly over all the [5 
traverse angles (Table III). 

2. Linear Errors of Closure 

When all angles have been corrected, the process of calculating the improved 
coordinates of all traverse stations may be done. The check on angular closure for a 
closed traverse does not guarantee that the entire survey is correct, because there can 
be considerable errors in the linear measurement of individual lines. Such errors may 
not show up in the angular check. In order to check the closure of the traverse, it is 
necessary to determine linear error. 

The linear error (LE), the departure error (6x), and latitude error (Oy) ina 


traverse are determined by equations 


LE = [ (8x)*+(dy)? ]!/2 (4.9) 


Ox CE eRGim (4.10) 
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dy = GN, - GN, (4.11) 
where 6x and éy are the traverse closure in departure and latitude, Giemeamd. Gris] ane 
the known and computed grid easting, and GN, and GN,’ are the known and 
computed grid northing at the closing station, respectively. By substituting the fixed 
and computed values of grid easting and northing in Equations 4.10 and 4.11 for the 
data of Tables II and VII 


$x = 325348.292 - 325348.180 = + 0.112m 
Sy = 3076179.297 - 3076178.924 = + 0.373 m 
LE = {(0.112)*+(0.373)7]!/2 = 0.389 m 


The relative error of closure provides a better assessment of the quality of a 
traverse than the linear error of closure. Therefore, it is common practice to calculate 
the relative error of closure, which is the linear error of closure divided by total 
distances of traverse, and to express the result in the form of a ratio with unity as the 
numerator. For the data of Tables II and VII, this computation 1s 0.389 / 14853.800 
or | * 35,185. 

Using the Compass or Bowditch rule, the computed traverse closures and 
corrections were obtained (Tables VIII and IX) and then the adjusted station 


coordinates (Table X) were computed. 


E. LEAST SQUARES ADJUSTMENT BY OBSERVATION EQUATIONS 
The adjustment by observation equations, shown in general form by Equation 

2.13, can be done more directly than the adjustment by condition equations. To 

achieve this, the explicit expression for the residual vector V from Equation 2.13 1s 


substituted in Equation 2.7 to obtain the following equation: 


(F 2 BX)! ew (F - BX) (4.12) 
= (F! - BI x! \wr - WBx) 

xT BTwex - x! Bl we + FIwe - Fl wsx 

x! pl wex -2F! wax + FIWE 


@ 


where X! Bl wer = FIWBxX are scalar quantities. 
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TABLEW i 
TRAVERSE CLOSURE 


I) Angular error 


Known azimucneat lace ices 293 5 Tieleeios 

Computed azimuth 273° 56' 49.58" 

Angular error | 22 

Angular correction per station Le 
II) Linear error before adjusting azimuths 


Known Coordinate ew) ae 
at Turtle 2 3253484292 SOy Oly shez 7. 


Computed Coordinate 

at Turtle 2 325348.180 
Error ox = + 0.112 
Linear error of closure 
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TABLE IX 
LATITUDE AND DEPARTURE CORRECTIONS 


Traverse lines X-correction (m) Y=correction( Mm) 
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The minimization of Equation 4.12 can be done by taking partial derivative with 


respected to each Unknown variable (X) 

®’ = 6@/8X = 2X!B! ws - 2F! WB = 0 (4.13) 
Transposing Equation 4.13 and rearranging yields 

(B' WB)X = BI WF (4.14) 
Or 

NX + U = 0 (4.15) 
where N = (B! WB) is the normal matrix of dimension u by u, and U is BI WF. Then 

Xeen lv (4.16) 


For the adjustment in traverse, the vector V in equation 2.13 is represented the 
residual of observed angles and distances. If there are n observed angles in a traverse, 
there will be n - 1 observed distances and the number of residuals becomes 2n+ I 
which includes the residual of angles and distances. 

There are two types of condition equations in the adjustment of observation 
equations: the angle condition and distance condition. 


From Figure 4.2 the angle condition can be expressed by 
Via = By - (45; - Gy) (4.17) 


where v;,, is a residual of observed angle, ai; is a forward azimuth, aj, is a backward 


azimuth and B. is an observed angle. Equation 4.17 is suitable when aii 2 a,. If ai; 


< d;,, the quantities in parenthesis must be added by 360°. 


The Equation 4.17 can be expressed in coordinates of the two stations as 





ea Xe (4.18) 
pie -] -1/2*k i 

¥; a heal) - ta 

la Bi oe! a a 
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Figure 4.2 Determination of Angle and Distance from Coordinates. 


and 


, 


me Sly aged ko 
gc [ tan oe Ae tial 


k i 


(4.19) 


When Imnearized by using Equation 2.12, Equation 4.19 yields 


Bi cam ce. ie Ni Yu Xp ae ole a OX; ar a50Y; i az0X; fa ayo; ale 


mere x. , Wee Np Yy. Xp and ae are the estimated station coordinates. By 


substituting the linearization of function g;, Equation 4.18 can be written as 
where 6X,, 6Y;,.....,6X,, and d6Y, are unknown parameters (correction in X and ¥ 
coordinates), coefliceints Apr AQseorrers ag, and a¢ are the partial derivative of function S 


with respected to X;, Yi... eee, respectively. 


lor the distance condition, it can be expressed as 


vig = dh 1(%- X.)* + ae y2 h/2 (4.22) 


l 


4] 


where 


hy = [(X - X.)* + (Y;- Yj y2 h/2 (4.23) 
The linerized form of Equation 4.23 1s then given as 


hy = hj(Xj’, Yj’, Xj, Yj, Xq’s Yq’) + bjSX; + bbY, + bzSX; + bySY; (4.24) 


Substituting the linearized form from Equation 4.24, Equation 4.22 becomes 
Veg) OO xeuaie > Olt gat b30X; -- b4oY; = inp (4.25) 


where Fl and F2 represent the constant form for Equations 4.21 and 4.25. Table XI 
shows the partial derivatives for Equations 4.21 and 4.25. 
Sometimes, the station positions which is determined by intersection from 
traverse Station 1s also to be adjusted simultaneously with traverse station positions. 
In this case only the angle conditions are added to observation equations of traverse. 
In the data adjusted under this thesis, there are 15 observed angles and 14 
oberved distances which consisted of 13 unknown points, consequently, there will be 29 
observation equations included 26 parameters of 6X and oY. Thus, they can be 


expressed in matrix form of Equation 2.13 as 


Vo9,1 + B29 26 X26,1 = F291 (4.26) 


where Vip V2 pore V15,1 are the residuals of angles; Yi6,b YIT,b Vo9 1 are the 
residuals of distances; B is the coefficeint matrix (consistings of a’s and b’s) of 
parameter X (Table XI); and F is the constant vector. 

When three intersection points with 10 observed angles were added (Figure 3.1) 
to adjust the coordinates, the observation equations have 39 equations including 32 


unknown parameters. 


Y39.1 7 2303p cpa econ (4.27) 
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TABLE! 
THE COEFFICIENTS OF ANGLE AND DISTANCE CONDITIONS 
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By using Equation 4.16, N7! is the cofactor matrix and the diagonal terms of this 
matrix gives the variances of the adjusted coordinates. To obtain the residuals of all 
observed quantities, the reverse process must be done. The correction of X and Y 
adjusted coordinates and their standard deviation (o’s) are shown in Table XII. And 
the adjusted standard deviation of observed quantities were obtained by multiplying 
standard deviation of unit weight (viwy / ry t/2) to a squares root of diagonal 
element of B(B! WB)! B! matrix (Table XIID). 
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ABLE 
ADJUSTED STANDARD DEVIATION OF ANGLES AND DISTANCES 


Number Angles Distances 
seconds meter 
+ ) ee 
1 .8 0.029 
zZ 8.9 OOS 5 
3 9.1 0.045 
+ eZ 0.024 
5 9.3 0.048 
6 9.4 O2056 
a 8.4 Ooo 
8 Sa 0.029 
9g 7.8 Omes. 
a) Sieg 0.034 
eae 9.6 OG. O0Z6 
nee Os), 0.048 
io 9.8 0.054 
14 ee O70 5G 
5 6.8 
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V. ANALYSIS OF RESULTS 


A. COMPARISON OF ADJUSTED COORDINATES 

If the given coordinates of the control points are assumed to be error free, then 
the accuracy of the traverse station coordinates depends only on the accuracy of 
distance and angle measurements. The adjusted traverse coordinates obtained by the 
approximation method are of a lower order of accuracy as only the errors in the 
misclosure in azimuths and distances were determined. These errors were distributed 
by assuming that all observed quantities had an equal probable occurrence. 

The least squares adjustment method provides a better approximation of the true 
value. Therefore, the adjusted traverse coordinates obtained by this technique provided 
better estimates for position of all traverse stations and the accuracy of the adjustment 
can be checked and statistically tested. 

After the 13 adjusted traverse station (from stations 2 to 14) coordinates 
obtained by the approximation method and by the least squares method were 
compared, the difference in coordinates at each station were computed and plotted 
(Table XIV and Figure 5.1). The largest difference was at station 14. Because stations 
| and 15 are held fixed, the least squares techniques adjusts simultaneously errors in 
azimuths and distances, while the adjustment by approximation adjusts errors in 
azimuths and distances sequentially. Consequently, the largest difference in traverse 
distances occurs at the last station (station 14) before closing of traverse at the fixed 
station. 

When the standard deviation of observed quantities before the adjustment 
(Tables II] and IV) were compared with those obtained through adjustment (Table 
ATI), the standard deviation of all observed quantities in Table XIII showed 
increments. That means, the estimated standard deviations were optimistic. 

In this thesis, the three intersection points were also adjusted. The adjusted 
coordinates of these were compared to NOS results. The largest difference occurs at 
point no. 3 (0x = + 0.140 m; dy = - 0.158 m). The standard deviation of adjusted 
coordinates at this point are Ox, = + 1.87, oy, = + 0.96 m (Table XV). 
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B. ANALYSIS OF THE REFERENCE VARIANCE OF UNIT WEIGHT 

The weight matrix was set for the least squares adjustment by using Equation 
2.11. The Te (a priori reference variance of unit weight) was assumed as 1. The 
result of Gy" = (viwv / r) was obtained after adjustment. The value of Gy* (a 
posteriori reference variance of unit weight) can be used to evaluate the weighting 
scheme used in the least squares adjustment. 

The standard deviation of the observed angles (6) and distances (64) used in the 
least squares adjustments (solutions |, 2, and 3) and the corresponding a posteriori Cr 
obtained for these solutions are listed in Table XVI. As the a posteriori ony forthe 
solution no. 3 is closest to the assumed a priori So" (= 1), the weight (or the standard 
deviations) used for observed angles and distances in this case seem the most realistic. 
Further statistical testing for cre done by Chi-squares or F-test was not carried out 


under this thesis. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

By using the weight of observed quantities for the adjustment, the traverse 
station coordinates computed and adjusted by the least squares observation equations 
method were more accurate than those obtained by the approximation methods. Even 
though the observation equations method may require a greater number of equations 
than the condition equations method, the processing of the data for adjustment 1s 
easier and the corrections in X and Y coordinates can be directly obtained through 
iterative solution. This method is suitable when a computer with a memory capacity of 
over 500 K bytes 1s available. However, for local work or a relatively short traverse, an 
approximation method 1s commonly utilized when economic and logistic criteria are 


considered. 


B. RECOMMENDATION 
The INDTRA Fortran program written for this thesis is automated for handling 


rs 


only two kinds of survey techniques: traversing and intersection. With the computers 
available at NPS, the development of adjustment programs for covering a wide range 
of survey techniques should be done to use and continue analysis of the mixed kind of 
survey techniques including traverse, triangulation, trilateration, resection, and 


t 


intersection. 
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APPENDIXN A 
LINEARIZATIONS 


This 1s a linearized form of m functions in n unknowns. 


Yy = Fy (Xp, Xgney XQ) 
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oh 6x 
Ym Pm 8m 
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Ax, 
AX = 
Ax 


The general form of linearized functions becomes 


_— 3a 
Y = 9+ JAX 
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APPENDIX B 
TRAVADJ FORTRAN PROGRAM 


This program is used for computing and adjusting the traverse station position by 


approximation method (Compass rule). 
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DETERMINATION OF THE GRID DISTANCES 
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WRITEC9, 180 
WRITE(9, 170 
WRITE(9 , 200 
WRITE(9 210 
WRITE 9 , 200 
WRITEC9, 1/0 
WRITE( 9,200 
DO 20 J=1,N-1 
MOIS = MDIST(J){OUMI 
RRO = ROC) }/DUMI 
S1GD = (ROC OMe ae 0100 
WRITE(Q 220 MNAME( LJ) ,L=1,2),(NAME(L,J+1),L=1,2),MDIS, 
*  MDIST(J),RD(J) , RRD 
CONTINUE 
MOIS = MDIST ON). OUMI 
RRD_ = ue DUM1 
MA TEC aes (NAME(L,N),L=1,2),(NAME3(L),L=1,2),MDIS ,MDIST(N) 
Siar 79 ( BRAN) *0. 00001D0)+0. 01D0 
WRTTECQ>¢0 
WRITE(9,170 
EORUAT CSE reaes 
FORMAT Ze KIKI KIA KKK KKK KKK KKK KKK KKK KKK KK KKK KK KK KK KKK KKK KK KK KKK KKK 
FORMAT(9K,'* : ‘ 
FORMAT( 9K, '* MEASURED HORIZONTAL ANGLE : ’ 
FORMAT(9K, I oko OK Ok kk kk a a ke tk DEGREES * MINUTES * SEC 
FORMAT( 9K, ie AT TO ‘ " 
FORMAT( 9K, _ ; : 
OMA pz 2A4,3X,"*", 3X, 2A4 3X tl SX ges one eamne ashe 
FORMAT( 5X, ###eRt LADAHBAHHHHAAAADEEHBHEPRHABHHRHHARARPHHRAARAHB HHH 
“Heda RHER HAH baat ) 
FORMAT(SX, # yy # # 
FORMAT 8X # 1 STANCES # MEASURED DISTANCES # REDUCE 
*) DISTANCES  #") 
FORMAT(SK, #4. # 4 # # 
FORMAT (SX 4 Fon # TO # METERS # FEET # METERS 


EORMAT( SY, ely 
*FQ 3 1X 
FORMAT 


gat #' 


2A4,'#',2X,F8.3,2X,'#',2X,F8.2,2X,'#', 2X, 
Ree) 


FORMATE OX. ee ac 
9 
* 


FORMAT(9X, '* 
FORMAT( 9X, '* 
FORMAT( 9X, '* 
FORMAT( 9X, '* 
FORMAT( 9X ,'* 


r BHAT, 3K. on 
FORMAT 
RETUR 


KNOWN POSITIONS 


2X,F 
OXi'*! "ay 3a Ax tx 


' , 5A4 


NAME 
OF 


i*l 
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* 


* 


COORDINATES 


x4 


* | 


FORO OO ok kk ek tL 


* 


* 


GRID NORTH ™* 


2X,F11. 3 
Loy? ety 3 2k eae 


1 x | 


* 


ix! 


eT 3s 


* I 


GRID EAST *' 


dey a 


C) ©) ©) CII IGE ICI IMO) 


CO 


C2 “G “Gc do )0 1010101090) 


CO 


) ©) Gr Cam C{eC 00010) 


END 


Sees eee OCC e ee 


@ CONVERT DEGREES, MINUTES, 6 
@ AND SECONDS TO RADIANS @G 


@ 
CECCEEEEEEEEECEEEEEEEECEECEEEEEEEECE 


SUBROUTINE DMSR ( D1,M1,S1,R1 ) 
REAL*8 OD1,M1,$1,R1,PI1,N60,N3600,N180 
DATA N60/60. 0D0/ ,N3600/3600. OD0/,N180/180. 0D0/ 


ODO*DATAN( 1. ODO) 
(D1+(M1/N60)+(S1/N3600) )*P11)/N180 


ee aga COCCI 
@ COMPUTE THE COORDINATES @ 
@ FROM ONE KNOWN POINT WITH _ @ 
f MEASURED DISTANCE AND ANGLE 
CEEEEECEEEEEEEEEEEECEEEEEECEEEEEEEEE 
SUBROUTINE UTM (X1,Y1,DIS12,ANG12R,X2,Y2 ) 


Reem cueAL, Yl X25 %2,01S12 ,ANGIZR,.DIFX,DIFY 


Die =" DIS 1i2* Beds ANGISR 
DIFY = DIS12* BeOst ANGLER) 
X2 = X1+DIF 
va = T1tDTrr 
RETURN 

D 


CEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE 


@ 
C 
@ CONVERT ANGLE FROM RADIAN @ 
& TO DEGREE,MINUTES,SECONDS @ 
C 


apeoeoeoeororoeorororererereeecs 

SUBROUTINE RDMS ( RR1,D2,M2,S2 ) 

REAL*8 RR1,D2,M2,S2,PI1,N60,N180, TDEG,DIF, TMIN 
DATA N60/60. OD0/,N180/180. 0D0/ 
OD0*DATAN( 1, 000) 


PIl = 

TDEG= (iBorse1)/ 

D2_ = DFLOAT vet NT(TDEG) ) 
DIE = wiWwEG=) 

TMIN= DIF*N60 

M2. = aa eS IND) 
DIF = TMIN-M 

S2 = DIF*N60 

RETURN 

END 
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APPENDIEG@G 
INDTRA FORTRAN PROGRAM 


INDTRA Fortran program 1s used for computing and adjusting traverse station 


position by least squares method of observation equations (indirect observation) and 


simultaneously adjust the intersection points which are observed at each traverse 


station. 


ne ee ee ae a 

ce . « a) & 

ce Path — Orr ec 

ea -— Wm & jz eC, 

ec WY OQ -—Lu ©, 

© me ee Mm Or ©, 

ce . OO © Lit Lite e, 

2, = OO eH -- o> e 

ce <l i 6G <t Li ec 

e, Gre Wm CO =z FUL 2 

o = Lu mm Zl e 

ce ~“ WY OQ <<oO ce 

e, = WY ce ox FO eC, 

eC, —) Lid Lu O NM ©, 
CCCCCe ec < ee oe Qa aL ec 
CCCCCEe ce Lit = 2 : Sac!) ec 
a ee ce = — ° Or ec 
ec ea ce -— fr CM RMlLUY ec 
a ae ce TreS= UWS ODTVZOO ce 
ec ea ec “—~ Ze Oto>LeEOZzzZ eG 
EC, Oo & C, F—- LUrRHO « OW ree EC, 
ec Oo @& e [p= a2Z—IOawy FUuLiL @OOrF-E-E-M CG 
eZ ur G& RD MQ reoQa WwW NHNOOxOZzTtDDnDeré CG 
Ge OF @ ae lu CO OO ZWwWIM<t Oo mmaakl eG 
e, Lt CG qe] AS 2 aeszsrtuwwdZzsm —O=Ss-t- G& 
aS F= G& me OO LOZ BODO OOZZHWOOW BG 
QoS = ae Qo Zr WS O@WZOrHSsD OCOO=z G& 
eS WM @& Qe ck NO] 2 Sota Oa z <= ©G 
SY =z CG eg (ee) WE perce le ps eettt=z oorcerr’ CG 
oe FO CGC eC, LWIMOOOAs«MiU LLaZzny RRS .<CHOOt CG 
CQ, NH QB NMHLWO ALLOMHALS>>za -FUUA FG 
oe Dre CG e& MAW +UWILIO Or WLI <tl<f eC, 
oe Dt & eS OZ+ AY BW COO CsarowsZz eG 
Gg AD» CG ec < — FeO Ui-Hrasz S  CUG 
eZ AVY @& © tUtLl--FZIULMNYORZOOO FWRHO & 
ec Uu & DB OOMMFVUOTOZY Cer lOMwFEZ 
Ce WwW eG 2, PHN Wetoodtatr OOD» ~~ & 
oe faz CG ee Ceroo rece SOOWOYOOO=Zz CG 
Co «OO BG @e, Lut WOOO 'CZZOrew erernm CG 
CQ, DH CG Ge OM. Se Ort Ota rnm> wom Q, 
ae Or fe QD SEO AODS CWE =F FrDOMMDDIL CG 
eye nt & CQ DBDDaoOwOea2aDVOe¢) || Maza AMD G 
ec —- & DR ZZOOrwNZ — O Il a 
eS Fe & ec run nooitl uid Wout & 
eZ MU & ~~ iiiuead il & icy ae 
ST an fc& &, Omnm aL HY == © 
eS d6 UM CG ae << SH OH xO OH ULL ©, 
eZ ADO @& Co 2a Or oe Ss OOo) —! e 
C, ©, C2, CTOOrFO NOOO <COFFAW2 <cl 2, 
ce C © 22222” 2142 Oe 22 (jc e, 
eo ec ae oe 
CCCCCSe e ec 
CCCCCEe CCCCCCCCCCCCCECCCCCCCCCeCe 


QOOdODOODOOOOOOOODOOOOOOOOOOOOOOOOOOOIOIS 


DOUBLE PRECISION 


x xx* x kK Xx 


DOUBLE PRECISION 


* 


x 


DOUBLE PRECISION 


* 


x 
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we CO 
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ae ” mist CO 
=ONIH OQ OIeo Qa. 
“<t <TH cr MIM 
Cw « <{ <<t Fl 
t-C> — = awn Gy 
Mf Se © eS —te> t- 
Mors « =) Li My <{ 
O-<LO SO AQ — 
“ aALie © Ce a WY 
Qs a Om > ce 
Nos - ON <t COD Lu 
o~ Cc ei-, =< ees Y 
SS IO NO ae ee oc 
co —>=<tt ~S oO Oa Lu 
Se nf ri >< — ——— << 
oO —LLONI > 4 <t OO ce CO 
© ZOn= =) ce — rH — ~~ © 
= “NI + OO -- Y) oo Ce we 
“ O--~~«<t_H a <C C505 UL. © Y 
o—~ ZOuUIc AI a C5 LULL oO LO YY) 
« OM a" cy) “ Lis 
© ft Zz — on 
CO OLN + -. Li 4 —~—™ Lu = Lu cr oO L 
Y~ 2000 COO oO }— NOY }— ” t- << re K =a 
>< =NIDt Oo ™O ce —_ <C Lu <= . = VW) = fiw 
Oo Cutt ONO WY <t QQ Zz i on © = ” =a 
O Ose ed JF 4D fo) = — ~ ao” & 
>< = ONIN OzO m4 N= oo CO = QO. — = = = 
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ses & #= s= 85 Se = a SS O86) Seicuaes 2705) 
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S56 4 - a eu = - Sone = o> Sas ae 
ee ct <O ~~ —~ =< Ui. Le. en 
OC ae Coe CI x x oO oO Ge crea = Zio "I eu on ae, 
oS - He <TC) a a MIM or <a rd 1Ocoor! — Smsn <= ui 
Ons MN LW CO Ss LoL <<t oe ee ee UI “MOMND=SWY 6 Or 
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READ(5,*) NZ1(1) 


. 
f 
: IF (RANG. EQ. 0.0) THEN 
C READ APPROX. GRID COORDINATES OF RESECTION STATION 
READE 2a} NCO 
: READ(5,*) ECO 
CNGRID(I) = NCO 
CEGRID(I) = ECO 
GO T10na@ 
; END IF 
K = K+] 
DIS(K) 7 Brees 
: STDD(K) = (DIS(R)*DUM11) + DUM12 
AZIMUT = AZFIR + OR they 
IF (AZIMUT. GT.A AN360 
AZIMUT = AZ Ni 
s END IF 
if COAZR(K) = AZIMUT 
Ceceee FOR PRINTING 
CALL RDMS (AZ MUT ,DUM1D ,DUM1M, DUM1S) 
CFAZD(K) = aN ea 
CFAZM(K) = IDINT(DUMIM 
CFAZS(K) = DUMIS 
ceceee : 
: CALL UTM (XGD,YGD,RANG,AZIMUT , DUMMY1 , DUMMY2) 
BEYER = DUMMY1-XGD 
; DELY(K} = DUMMY2-YGD 
SUMDX = SUMDX + DEL 
: SUMDY = SUMDY + DELY(K 
RGR IBC T= DUMMY 1 
, CNGRID(I)= DUMMY2 
C 
: CALL GRIDAZ (DUMMY1,DUMMY2,XGD,YGD,AZC ,DUM1S) 
C@@CRE FOR PRINTING 
: CALL RDMS (AZC,DUM1D,DUM1M,DUM1S) 
CBAZD(K) = cae: 
CBAZM(K) = IDINT(DUMIM 
CBAZS(K) = DUMIS 
ceeeee 
AZFIR = AZC 
EO KE = DUMMY1 
: YCOY(K) = DUMMY2 
E 
XGD = DUMMY1 
YGD = QUMMY2 
, TODIS = TODIS + RANGE 
C 
C 
C 


100 CONTINUE 
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E K10 = K=1 
beet =] K10 
ane. . = feast} 
eC =" YCON! I 
“3 CONTINUE 
C 
. KK = K+l 
READ (5,105) D.M,$ 
ea any tras 
MW = DFLOAT(M 
i SS = DFLOAT(S 
READ (5 aa RANG 
: DIST(N) = RANG 
c 
; CALL DMSR (DD,MM,SS,R) 
eeu = 
; STDAC(N) = DUM1O 
AZLAS = AZFIR + ANGI(N) 
IF (AZLAS. GT. AN360R) THEN 
AZLAS = AZLAS - ‘AN360R 
; END IF 
DL RAzT Som A 7a! 
WPRED(1)= DIFAZI 
: CORAZ “= DIFAZI / OFLOAT(KK) 
COeeRE FOR PRINTING 
: CALL RDMS (DIFAZI,DUM1D,0UM1M,DUM1S) 
FA710 = IDINT( ABS eile 
FA71M = IDINTC ABS (DUMIM 
: FAZ1S = ABS (DUM1S) 
: CALL RDMS (CORAZ,DUM1D,0UM1M,DUM1S) 
FAZ2D = IDINT( ABS ae 
FA7Z2M = IDINT Abs DUM1M 
FAZ2S = ABS (DUM1S) 
ceceee 
WpREDGS) = SUMDX - ( FaRTOC3} - ee 
; WPRED(3) = SUMDY - (NGRID(3) - NGRID(2 
DIFDIS = DSQORT( (WPRED(2))**2 +. (WPRED(3))**2 
ERR = PBINTE ODIS ¢ Atepes | (3)) 


C 
c READ NUMBER OF MEASURED ANGLES AT EACH TRAVERSE STATION 


ee READ( 5. Ny N2_ 
Numéa) “= 


101 CONTINUE 


Oo 
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TTCOX(II) = DUMMY2 


¢ = 
ee CONTINUE 
C 
C 
; ASP = NSTA - 1 
C@@ KEPP APPROX. GRID COORDINATES OF EACH STATION IN 
wes ASNY AND ASEX 
DO 102 I=1,ASP 
AREY ET = EMER IDET 
ASEX(I) = CEGRID(I 


102 CONTINUE 
C@@ NUMBER OF OBSERVED ANGLES PLUS OBSERVED DISTANCES 


NT = N+K 
C@@ NUMBER OF OBSERVED ANGLES 
NA = N 
C@@ NUMBER OF OBSERVED DISTANCES 
ND = K 
C@@ NUMBER OF TRAVERSE STATIONS NOT INCLUDED KNOWN STATIONS 
NP = ASP 


C@@ NUMBER OF UNKNOWN DX AND DY 
NP2 = K10*2 + NR1*2 


“oO 
‘a 


NI = NP2**2 + 3*NP2 


NUMBER OF TRAVERSE_STAT 
PLUSE NUMBER OF INTERSE 


NC = K10+NR1 


@CEECC CALL _SUBRQUT 
BY LEAST SQU 


CALL INDLSQ 


* 
* 


Fone eee ee mg os CP ear t weet SEG 
“a 


o~ PR 


@@E@ PRINTING DETAILS 
TEE DEVIATION OF UNIT WEIGHT =',STDAD 
WRITE(6,*)'STANDARD DEVIATION OF ADJ. ANGLES! 
0 102 T=1,NA 
CALL ROMS CVC 1 DUM1D DUM IM , DUMLS) 
MM Ps DUMID*3600. ODO + DUMIM*60. 
TTE( 6. =e ADJ. STDA. 


CIM OIOY 


103 CONTINUE 
WRITE(6,*)'STANDARD DEVIATION OF ADJ DISTANCES' 
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=" VM(1,K4) 


et, 1 ADJ. SDD. 


105 FORMAT(I3,1X,12,1X,12) 


104 CONTI 


c 


i on™ HOO cS 
a =a So 

=—C) ZN At <Z 

oe CALICIZZ™ eS 

>a OFAN = 

OO > =O Zot Z 

OO -« pet tons 
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mone D OZ=zale4 oe 
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O an « OVE NEY me 

Keo 4 aia SS Ole 6 — 

SS > aw, S>O = 

HY => SZC DO) HH ™ 

Kc) an on Ow m0. NSF <Q = 

Naz Ss «OL «8 a Oe oO 

— O_- NS ZI Li} es 

Ont -a_J to 4—>.>-<f Saas . 

MALL An Se Oe © 

Snr O OD 4-2 > > Ra 
CCC0CéCCéE8, we ae oN Se Sz © 
CCCGZCCG MSM FQ WAH OF a= = 
©, @, Sie Om a on WZ HI —2 
©, @, “ial LU aN Ow = 
©, © NI SE MOA 1 “ 
@ 2 =<IOZWN Oa = ae Bie ~s 
©, Q Da DOD -- Sony aL 2 © 
2 @ Zeb = INTO roQ Cc) 
ec Cc Zn => O<ttua.+--O< O-s¢ © 
2 a mM ~GeS Sp ee a = © 
© e NANA O9 Myx LAWN oN S 
© wu @, Zama Oa 4-OQ *« Lid cy © 
CG, Mr Q, "ZA WS OOS yee = 
eo xO CO, Oo Oe <t “LO. 0O-™«L ae © 
qe Diu Qe ZA. Ss «a aA a Ola = ny = 
SG. Gc 2 "2 stan) >—-N_OAIONY AY A a 
C, We @, oa f. = = HOOK Z<te>- ae CO 
@, e) 2 OOF =M Sa o_O “a © 
Oe +z © a 2 O/C DDNWIS I SO "Ks OC 
Qo Ne 2 — OO<th-n>Fr- ze . 
Qe << ©, CO 
em twtr ©, Bg = 
eZ Hon CG wW Pi — 
©, = © -J (es) lam 
eo uo G OQ ==! ” 
So 2zze G&G a W” iu 
C eitUt- G — s rs 
eZ F-=et CG i 
eG D> @ Lu WW) = 
oS Once Ge Za oO oO toe 
Oo wal & — Qa. “NN 
eZ MDD” CG = ~~ a 
eo DOM G = Lu Lu 32 Heat a 
eo nO @ © 3 C5 
©, @, oc ao Lu << 
e, Q, aa) =) ae = SiN 
2 ©, = © ee <t Wa. 
CCC, ~ O = a) ae 
CCCCCO,Q, xxx * xX KRXKXXX* x 


OOWIOSSOCSMOOOMOO © CoO Oo WO 


69 


1 = STDA(T) 


THEN 
1,NA 


) 


PW. NE. 1 
O 40 I 
DUMMY 


IF (F 


C@@ SET WEIGHT MATRIX 


Ce@e@ FOR STD. ANGLE 


g 
C 
C 


CALL DMSR_ (NUMO,NUMO,DUMMY1,DUMMY2 
= CHBot Ts = 1 cna ad te ENCURMY2) )**2) 
40 CONTINUE 
C@@ FOR STD. DISTANCE 
Does oma 
J = NAtI 
DUMMY1 = STO0¢T) 
STOO(J)= 1.0D0 7 (DUMMY1**2) 
50 CONTINUE 
C0e SET UNEQUAL WEIGHT 
DO 70 I=1,NT 
ite TEL 60 J) THEN 
VIN P-3y 2 STDO(I) 
END I 
60 CONTINUE 
10 CONTINUE 
END IF 
C 


aes PRINT WEIGHT MATRIX 


IF (PRINTD. NE.0) THEN 
RITES, *) {WEIGHT MATRIX’ 


SWFM ('R-C. ',NC4,WM,NT,NT,NT,NP5S) 


99 CONTINUE 


K1 = Kitl 
[Eki oie g) THEN 
e223 
END IF 
CHECK = 0 
CALL SUBROUTINE TO CALCULATE "Fl" MATRIX 


MOM «=O 
@ 
‘@ 


CALL CALFM (NUMM,K 
* NT NGRID 


MOIT1IO0IOY 
@ 
@ 


CALL SUBROUTINE TO CALCULATE "A" MATRIX (COEFFICEINT DX,DY) 


CALL CALAM (NR,NZZ,NUMM,KK,K,DIST1,TCOX,TCOY ,NT,NA,ND,NP, 
* NP2 NGRID,EGRID,ASNG,ASEG, BM) 


3©GI63 
‘@ 
@ 


TRANSPOSE OF "A" MATRIX 
DO 120 I=1,NT 
DO 110 J=1,NP2 
BMiCdel,) =e eicleap 


C 
110 CONTINUE 
120 CONTINUE 


C@@ "A" TRANSPOSE * W 
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CALL VMULFF (BMT,WM,NP2,NT,NT,NP2,NT,BMTWM,NP2, IER) 
cee "A" TRANPOSE * W* A 
CALL VMULFF (BMTWM,BM,NP2,NT,NP2,NP2,NT,NM,NP2, IER) 


C@@ "A" TRANSPOSE * W * F 
: CALL VMULFF (BMTWM,FM,NP2,NT,N1,NP2,NT,TM,NP2, IER) 
C 
C@@ (INVERSE OF ("'A" TRANSPOSE W * A) ) 
: CALL LINV2F (NM,NP2,NP2,NMI,N1,WK10, IER) 
¢C 
C@@ CALCULATE DELTA MATRIX (COORECTION VECTOR) 
C  CINVERSE OF ("AY TRANSPOSE W* A) * VAN TRANSPOSE * W* F ) 
: CALL VMULEF (NMI, 1M,NP2,NP2,N1,NP2,NP2,DELTA,NP2, IER) 
C@@ UPDATE APPROX. VALUE 
¢ 
C KON = K-1+NR 
DO 130 T=1,NC1 
TOKE) = TOOK 1) + DELTA (E-1) 2 + 1) 
} TCOY(I) = TCOY(T) + DELTAC(I~-1)*2 + 2 
130 CONTINUE 
¢C 
II1 = (K-1)*2 
: DO 131 f=1°NP 
IF (NZZ(1), . EQ. 0), THEN 
Azeal = ASE (1) + DEL TAK (171}%3 + 1) 
ny ASNGUL) = ASNG(T) + DELTAC (1-1)*2 + 2 
DO 132 J =1,NR 
TF (N2Z(1) . EQ, J), THEN 
ASEG(T) ‘= ASE (s + DELTA it(dr}}3 +3} 
ASNG( 1) = ASNG(I) + DELTAC If1+(J-1)*2 +2 
GO TO 131 
END IF 
132 CONTINUE 
END IF 


C 
os CONTINUE 


itete 1020 ‘al 
WRITE(6, 1030 
CALL USWFM ('R-C.',NC4,FM,NT,NT,N1,NPS5) 
MRITECE: 1088} 

CALL USWFM ('R-C. ',NC4,BM,NT,NT,NP2,NP5) 
VRETE(® 1088) 

CALL USWFM ('R-C.',NC4,BMT,NP2,NP2,NT,NPS) 
VRITEC® 1800) 

CALL USWFM ('R-C. ',NC4,BMTWM,NP2,NP2,NT,NPS) 
VRHTEC® 1098) 

CALL USWFM ('R-C.',NC4,NM,NP2,NP2,NP2,NP5) 
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WAITEC8 1080) 
CALL USWEM ('R-C.',NC4,TM,NP2,NP2,N1,NP5) 
VRTTECE 1888) 
CALL USWFM ('R-C.',NC4,NMI,NP2,NP2,NP2,NP5) 
WAITEC® 1790} 
WRITE(6,1100 
CALL USWFM ('R-C.',NC4,DELTA,NP2,NP2,N1,NP5) 
WRIT erertops 
; WRITE(6;1110 
DO 140 I=1,NP 
WRITECe ) T,' N= ', ASNG(I),' E = ',ASEG(I) 
WRITE(6.* 
140 CONTINUE 
END IF 
é 
Coe CHECK DELTA MATRIX 
c 
DO 180 enue 
- TEST DABS (DELTA(I)) 
IF (_ TEST2 .GE. TEST1) THEN 
CHECK = 1 
END IF 
180 CONTINUE 
c 
IF (CHECK .NE. 0) THEN 
0 TO 99 
END IF 


c 
C@@ CALCULATE "A" * X MATRIX 

CALL VMULFF (BM,DELTA,NT,NP2,N1,NT,NP2,STDO,NT, IER) 
C@e@ CALCULATE "V" MATRIX 


DO 190 IIeNT 
STDO }) STDO tT} - FM(T) 


eae STDO 
190 CONTINU 


CORRECT OBSERVED ANGLES AND DISTANCES 
DO 139 I = 1,NA 


MOOD 


ANG(I) = ANG(I)+(STDO(I 
139 CONTINUE. (1) +¢ (1)) 
DO 149 I = 1,ND 
NND1 = NA+I 
DIS1(I) = DIS1(1)+STDO(NND1) 
149 CONTINUE 
IF (PRINTD NE.0) THEN 
RITE(6, 1000 
WRITE( 6.1120 
CALL USWFM ('R-C.',NC4,STDO,NT,NT,N1,NPS5) 
WRITE(6 , 1000 
WRITE(6. 1130 
CALL USWFM ('R-C.',NC4,VM,N1,N1,NT,NPS) 
END IF 
C 
c 


C@@ CALCULATE "V" TRANSPOSE * W * VY MATRIX 


ee 


TNT. FM,NT, IER 
NT SSTANO NIE TER) 


NT.NT,NI.N 
NT NING 


Nd 


’ 


WFM ('R-C.',NC4,FM,NT,NT,N1,NPS) 


DO 


VMULEE aie 
VMGLFF (VM, FM 


} 
? 
E 
V MATRIX',/) 


= STAND * (DSQRT(DABS (NMI(I,J)) ) ) 
NMI(I,J) 


J) 
R~C. ' NC4,WM,NT,NT,NT,NPS) 


= STAND * ( DSQRT( DABS(WM(I,J)) ) ) 
x 
A 
A 
R 
R 
A 
/) 


a "IH - SSFP ws eS xe 
eo Sa NUE iftaqh xD>>e> 
t= 


CINVERSE "A"TRANSPOSE*W*A) * "A" TRANSPOSE 


oo SS OC. ee P< DK De OC OK OK OK OK OC OK OK OS OK OK 
San = LOLOLOLOLOLOLOLOLOLOLOLOLOLO 


DSQRT ( STAND / (DFLOAT(NT-NP2)) ) 


N 

J 
(I, 
ele 
Tid 


END I 
UE 


= Sul ee 
—— EE FF EE 


CALL VMULFF (NMI,BMT ,NP2,NP2,NT,NP2,NP2,BMTWM,NP2, IER) 


CALL VMULFF (BM,BMTWM,NT,NP2,NT,NT,NP2,WM,NT, IER) 


STAND 

F 

N 

0 

N 
RETURN 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 
FORMA 


i 
8 
a9 
Q 
2 
3 
- 
5 
6 
i 
8 
g 
Q 
] 
2 
2 
: 
9 


C3. €) Coc) Co) (>) QDI MaMa aa amt et atest et et ie 


O 
C@@ CALCULATE "“A'* 


OOo OO Oo O 2 


iD 


CWI IO IC I IOI IMO) OI) O09). O90) Ok 


™ 9090 ©9009 ©2909 


bop 
“IOV 
OO 


FORMAT Oe 
FORMATC/// 8X ae TRVERSE A TRANSPOSE *W*A)*A TRANSPOSE! ,/) 
END 


ELCEELEEEEEEEEEEEECEEEEEEEEEEEEEEECEEEREE 


CEE SUBROUTINE FOR COMPUTING CEC 
GEE F MATRIX EEE 


EEE 
EECEEEEEEEEEEEEEEEEEEEEEEEEEEEEEECEEEEEE 


, SUBROUTINE CALFM STA 
i NO 


DOUBLE PRECISION coca G 


* 4 
Ora 


INTEGER T1,12,K,KK,KKK,N,STN,KOUNT( KK) 
DATA N360/360. OD0/,NUMO/0. ODO/ ,N/0/,STN/1/ 


CALL DMSR (N360,NUMO,NUMO,A360R) 
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SUBROUTINE GRIDAZ (X1,Y1,X2,Y2,ANGLR,DIST12) 


REAL*8 X1,X2,Y1,Y2,ANGLR,ANGO,ANG9O ,ANG180 , ANG270 
REAL*8 ANGOOR,AN180OR ,AN270R ,DIFX ,DIFY,DIST12 
DATA ANGQ/9. 00/ ANG90/90. OD0/ ,ANG180/180. 0D0/ 
DATA ANG270/270. 0D0/ 
CALL DMSR ( ANG9O,ANGO,ANGO,ANG9OR ) 
CALL DMSR ( ANG180,ANGO,ANGO,AN180R 
CALL DMSR ( ANG270;ANGO.ANGO,AN270R 
DIFX = X2-X1 
DIFY = Y2-Y1 
IF((DIFX. EQ. ANGO). AND. (DIFY. EQ. ANGO)) THEN 
ANGLR = ANGO 
ELSE IF(DIFX.E ANG) i 
IF(DIFY. GT. ANGO) THEN 
ANGLR = ANGO 
END IF 
ANGLR = ANI18OR 
ELSE 
IF(DIFX. TANGO) THEN 
= ANGLE = ANGQOR-DATAN(DIFY/DIFX) 
ANGLR = AN270R-DATAN(DIFY/DIFX) 
END IF 
END IF 
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SUBROUTINE DMSR ( D1,M1,$1,R1 ) 
REAL*8 D1,M1,S1,R1,PI1,N60,N3600,N180 
DATA N60/60. ODO/ ,N3600/3600. ODO/ ,N180/180. OD0/ 


(ji QDO*DATAN DO) 
(D1+(M1/N6 }+($17N3600))* PI1)/N180 
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SUBROUTINE RDMS ( RR1,D2,M2,S2 ) 
REAL*8 RR1,D2,M2,S2,PI1,N60,N180,TDEG,DIF,TMIN 


DATA N60/60. OD0/ ,N180/180. 0D0/ 


PIi = ODO*DATAN( 1, 0D0) 

TDEG= (Ni80*RR 

D2_ = DELOAT( ID T(TDEG)) 
DIF = TDEG-D 
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$2 = DIF*N6O 

RETURN 

END 


83 


LIST OF REBEIE NCES 


e =. 


Computations, Fifth Edition, Heinemann Publishing Co., 1968 
Bomford, G., Geodesy, Fourth Edition, Clarendon Press, Oxford, 1980. 


Brinker, Russell C., and Paul R. Wolf, Elementary Surveying, Sixth Edition, Harper & 
Row, Publishers, 1977. 


Clark, Late David, and J.E. Jackson, Plane and Geodetic one Volume II Higher 


Allan, Arthur L., J.R. Hollwev, and J.H.B. Maynes, Practical Field Surveying and 


Surveying, Sixth Edition, Constable & Company Ltd., 1976. 


Davis, Raymond E., Francis S. Foote, James M. Anderson, and Edward M. Mikhail, 
Surveying Theory and Practice, Sixth Edition, McGraw-Hill Company, 1981. 


Department of the Army Technical Manual, 7M 5-241-8, Universal Transverse 
Mercator Grid, Headquarters, Department of the Army, 1958. 


Evett, Jack B., Surveying, First Edition, John Wiley & Son, 1979. 


Federal Geodetic Control Committee, Standards and Specifications for Geodetic Control 
Networks, Rockville, Maryland, 1984. 


=e . . ° 


Measurements, First Edition, Van Norstrand Reinhold Co., 


Mikhail, Edward M., and_F. Ackermann, Observations and Least squares, First 
Edition, Dun-Donnelley Publisher, 1976. 


Schmidt, Milton O. and William H. Ravner, Fundamentals of Surveying, Second 


a, 


Edition, D. Van Nostrand Company, 1978. 


Mikhail, Edward M., and Gordon Gracie, Analysis and ete of Survey 


Torge, Wolfgang, Geodesy, First Edition, Walter de Gruyter Berlin, New York, 1980. 


84 


U2 


Cay 


INE IAL DISTRIBU ITOMW ES 


Defense Technical Information Center 
Cameron Station | 
Alexandria, Virginia 22304-6145 


Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93943-5002 


Chairman, Department of Oceanography 
Code 68 

Naval Postgraduate School 

Monterey, California 93943-5000 


NOAA Liaison Officer 
Post Office Box 8688 
Nfonterey, California 93943-0688 


Dr. Muneendra Kumar 

Defense lapping Agencv. 
Hvdrographic/ Topographic Center 
Geodesy and Surveys Department 
6500 Brookes Lane 

Washinton, D.C. 20315 


@iiice ol the Mirector... . 

Naval Oceanography Division (OP-952) 
Department of the Navy 

Washington, D.C. 20350 


Commander 
Naval Oceanography Command 
WOE Mississippi 39529 


Commanding Officer 

Naval Oceanographic Office 
Bav St. Louis” 

Pode ississipp! 39522 


Commanding Officer ia 
Naval Ocean Research and Development Activity 
Baeoe, eouls- 

NSTL, Mississippi 39522 


Chief of Naval Research 
S0OUN? Quincy Street 
Arlington, Virginia 22217 


Chairman, Oceanography Department 
L.S. Naval Academy 
Annapolis, Maryland 21402 


Topographic Sciences Department 
Cate ESDENIC 
Defense Mapping School 

Ft. Belvoir, Virginia 22060-5828 


85 


No. Copies 
2 


13. 


14. 


ey 


CAPT Glen R. Schaefer 

Code 68 Se 

Naval Postgraduate School 
Monterey, California 93943-5000 


Supote Klangvichit 
Roy thet Navy Department 


Den ekoka at Bencrer 
Thailand 


Hydrographic Department Library 
Hydrographic Department 


Roval Thai Navy 
Bangkokyai, Bangkok 
Thailand 


a-Te 
ac 1343. NPS 
ene ae ‘California 93943 


Wang, Chih-P1 we 
SMC 1342, NP 
RiGniCiae ‘California 93943 


86 


17898 








vi 


_—— 
a 


WAT As PO 
HovTEp. 


Liog LIBRARY 
STGP ADIN TH 4 SCHOOL 


uY, CALIFORNIA 95946. “3002 

























5 Lah ed 






idee SM ae ee) oe eed ee eS peri t Prorty toe yi at 





thesK57282 Ba 5 Ven - 





ee ergy Fara Robs apes oom nd A Ere Ayre it ty te : 
Nie rs pete RES 4 CL a a cae ee eT eer Seay e. oy ey aries) , 5 7 : i . 
} nd a z JRO BMD. Cob Sagi see Gh pres Fy A 4 
(one ah rere oP ts Over) we Cura te ee ee ty Py oe Ard Baik dan cetr i TUT ats r re ; i 
spn fd) ede ce Rg Sa ee eee ave Ae oe ay oe oe ee eee ee , e] : t F . a 
are) 


“SAid bk SA et peued of MeaS Se A ey BL eet COTS Cee Pete te Py oe ar) 
Ft pte oad are Ze Ret eae Were OO cet ee er ee ee 


i | i WHT | 
p yt ioe FED rsrhOe EAA WN A Teh DDE DTH D WIR hd PDS heh Ate As (il Hi les ae | 
Ce yom Lo be Po See ee cr ae aCe Pe Pee eee: aw | mur Ue 5 
Pde oe Leeks re err es pear ae re eet) a ye Rae - eer } A rt 
POT Oe re ee ee sete. eS TS ee ee te | . were oo ue ary Hi y Fi ae fi 
METH 



































































































































































































Lbs ae det tee ee A er ee ee em re a tee ar Dis sh rr ed® Ba ** rl at 
ERE 8 0 shh n2 ED A mma hDREUL DS 4 18 +) el ee ey Tie we eee i ie ee Te ® 
Pbitetd Wi Seine CB ee LT RP Pe eee Reb oe Orrbeio.e | athe -vod ot 1 e bi 
Shao) Die LOL eS tT ae ae retire eerily er Male ee wa | were | i ' oT ‘ 
Gta adk Ole eye han ow Waeten ) a me 1) h s a et a ys A 
Lal seta Ramee 0.8, Tee er Lee Te ur cea "7 b ry Pr P r r P 
edhe eal Meet ad Be d Ppt wir ee an) i ? f 2 
aks Ble The PAPA kt akon 6.5 Per") rr . Fi rn , 4 Py 
Et TE Ott Oe eee et errs pred | ae: a P P ry Fi “ 
rt ee ty te me oy CoH .0:8 OAS mel vechirt Aehr-k cib yee xX LIBRARY 
ee Lt) Se ee Ce ny Ly | - PP oi eels ee reed T DUDL - : Py My 
err har te Le ate a Tice Pel oer A = r P P 
Pray 2) + cy eee be € Pe ay Cree P 3 J : i 
it We aed Ste eT 4 ea hyd ee phil Pip ba ~~ 3 Py ps ¢e.f A . Pa) ’ 
oh Tae he Ld A re r ee Lee Paw nt he ty ee a oe Yee eee ie NY er Teer eee | at ad | 2 f r) Ps ® P| F 
Intense Pe ty oa een ae Pht Prati ee Ly Lea pee Peery y an Pye te Ti ae oe | . a Oe S oan) , get f a) 
rer ieee SY ie. re ye ve ee Pe re Rawat ooh Pele Bee nar , cP ry sah 800:4-046 6 ° * E ~ os > r - Fi 
Ce oa Per cr entity a Seay a te ey Pr i in foe re teh: iat Pate a, é S ‘ ch ¢ ; 7 r 
3 ested dan sa xe ph AEA 1 Ae Ue ee Ye ty el a ef ry $b Pad ed  EBrckr bed Es + , at ere te | s an " 7 , ; F 
pele Eh Late a Ce ot Say a Ott oe Seley eran 7} Perry iri yee ee Aaa we irr yt er wie any. ree P t Part 17 28 , Fi Fi 
ame ld doe Me Mie ae YY eee ete ey! ar ee ary ery ee hee O sects Bete kus ¢ at ci ed ] s te ft | , ry 
os oe ee a ee Pe) O96 Reem ban Obekne Botnet 2° re cass Any bar a ee ae ere reer yee) aa Ey) P A 1 P 4 P} Py 
bP ier ae ghelae SOT bPAT RT Eke Ae ee é ee ear Sreh peed rae es TST bee a . < 8 4s ‘ ‘ Ca " 
gph tebe oT es ee eo a: ot 4 Ge Ps p 8 degde eri 4 pee aT , Bidmand £ of fF c 
eT ee ati ay ead “> u oe sod Ae Ye. ee ay eee x TTT ery eo) er o@er ‘ Pee 4 ere i) p ry » P 
cds sa aban nia c TEL ot SU aes oe Oe Co ae ee ee re ty ce Yor fey rae eee a ee ie i & 4:6 A@s a | s f ' 
daar he cid etd tata bed Pewee rot yar ee Pe. Pat) err oe ae ee | aud thee ay 7.2 wed . F a er ee a? oe c F A F 
ee ee et a eT Crs To On ey oe ay ey Pe Lea rae wr Oy Pe mee oe Ls 4 PST PL be ae ee in Py rae rr f Pi Fi , Pn Py ° 
ede eee ee ae Aimfaage® ae Perriit. rear: ret a On at re ee ke a Pee ere eee Be Py oa a ee D ’ , " | an) , ‘ H 
Palate ay Pr bs Ord ghee torte, r) Te a ee a ee Oe rer we rar) Po) ei 4 i ‘ ad Py PT r Py ' r , " 
en Ye | am Pe ee ee Se ee oe ee. 2 ar a Tvend AOE MIRA MOR GBs Fs © 8 Mey INR d ont > ' oT ‘ Ce ee ee pe ae | Py A , 
ret eka) a oe SY oe oY Be Peru Sar ® Bal, ae (oer oe tee Mare @-s°e - rr ae i ee | ’ td a ’ i ” e TY ri) ‘ 
i oT "ne as se eshass 44 Ps Pav A Tae} an ay Cao ee ee oe eee sof Sa it ae at Ci | D Por # ae | P A 
it ad dd a) ey et ee | a ee (ewe reer ee er ery Mert Cie oleae nee) 4 Pee ee ae n 3 PPT er ieee aera ey ee ‘ ‘ Pa ae | ry fs . Py ‘ e F ‘ r 
torte ein ESN ? Bons Mek i lgee fhe Siege pity jifyb ine 6 vee Rt a? eee ee ay) ’ er ae Se Y ; ' F , " p 7 . 
Pf 6 PoP. Bn @,V00 " ; FO AERO PCE TE een ee deeae it Paea@eetigagderan er r , F : 
ee Utd See! Tei er rey See Re Le Oe ee rere oe aes Pi ee ae We ae , $ Me Tarte. s * 8 ts ae ‘48 ae 
Tas ae ee ee Bt Puma } Setnent ons Perm y, a ee er | [ 3 a fo ’ L 
ba : Pot ne der ee et eee ee ee Pr eee Fr : err Py 1 © * per ee ‘ = 4 Paar Pa ‘ 
OE Le ee Pie ee ee a ts ee ee ae ee fod Bie Rd ALS Shee PPL een: ; ae at wor ry sus are F ’ r ~) F 5 
ee Se eS Oe onl oe Sa ee ene be ieee eee ae er eT ee avr: Py are F Fi ‘ f Ss = 
ly AO a ene et eee IT | Te ee Anka ct oa8 ip Pr Sas AB phon der te ofa ss zo45 , ee 7 a | i ee | ry Pi Fl 
Wy LAE Bed tg sos rrr ae oe for OF Bm HS te ra mer Fa iC b ‘ OPt J ’ A = 4 Py 
Fe Rogie | CANS R AFR OR Om f a rir Ca & * " 8,a 2,4 - ry | ‘ ® is Pr r n ; o ' 
Te a ae” rt ee er ee, | | Per) a ie ei ee Pe are UT ree 9 er Viper nt Eo 3. 0 We ca a Chiat] o , 4 i s 
hd ated dd 1 rae) Cee er ee | oa) A 1p gifey & ye ee ee ; Pay ay eo ar - 4 ee eae eee 4 r 5 A , = F 
‘mek wool ead t ON OA Toe a eee) Thee ee et i er os 1 D ro Pr ae a Ty F rr ‘ ae ae Ca Tae ae Fi a rr P ¢ fn ‘ ‘ D ‘ P ar Cr 
€ ROR R Semnh vo, waderety met s rn oe oer er ae ae rae ree erie ee CY ee oe ee ae ee ed O eae r en 7 r ag 1 é La ' i hs 4 
ad Nod J Te eee es Or ae rat ae 2 a ry on ’ ry 7 ' ret ° Ss Py F s 
eEaG tee fine arnt Aw oAneehy eRos.t os Prey oe tay Ag OAS eer ee tT ee a PTT Rr aaa were FY aa r ae a ae i ” r] ’ A 1 5 
tee Ce oo eae er ee a | oat pee A ee Pe) eve oe ee ee ee ae rd "{oarpio ae Btegpe Pa ae | ry id t . r : 
Ce a Pei ere tee ee St ee es ee eye oe er ey we ee ei ek ed rg Rieeg f Udtud Sob "a gud F P ri r 2 ,?, 8s Py 1 . ‘ ‘ 5 
Ud ee Ce ee Soni dyorm _¢ uit Cr On wre ae ae eee Re Ly eee ey r 7 ae Fo ee Te r ' un « L tt 5 " a ' bs 8 
oe ) ee cee 2) wre a wy ia) eres a Lie Sa eee oe ee © Obes R86 OA Uap ekee * nar aif tite : Pore 4 F F F F 
Ler Lat ae ee Pe 4 MEX. Ah re! ood PPL ese ek we ee ey aT ee 19 ee re Sr Pies | ry rear ier TY yy ie ere Sr ee , R 5 P - F eu rf 
hes a oS | ee) ee Pe eee) ee Se et ete ee ee ae PY Use Ce Pee ee Oe ey Pete eee ee A ted ital ey A n et tar s¢ " A F ene 
Ce ee er ee ee ee ee Po oe ee ee ee ee Dee ey ee Pt Pet ore rae Geen tt gar.e pie 5 t Cee LU ‘ a U iu Ff x = 7 
UE Ce ery Oe Pee re ee re Paver Ree Sree Pr ay Ye ee ee ee OY) Pye Be Cet ae er er ae PY eres ee de ee Le r) 4 ie o ai Fl 3 r a Hy e ’ ' My 
seals © Fd eae int ee ee eee ere) ay ey aa ar ae Ps a ee Or. Sie ee Of ey ‘ rr Aw ‘ ¥ p.cua ml F A A 
pe | Stet ah Haves ey ee er Pe ee Coe er, i 2 yt oe 2 ye ae ie ey es ere oe Te ee ee e | e ‘ Po ‘4 8 5 5 
Shad ne TC tt ae] ae ee ee Poe el Pea ere Aor ee Poe ee be Uae ie ee ra a ey ee ove PHRopu t . A P Sic ar as Pa ava _ 7 | = 
ae vo ee ae rar eer er ae se Ta oe a Pett? ee Pee) ae ee 41, a rafres 258 op idis caps as one Te en -pes ra ae : ri § r) fn ) 1 oe ar ’ 5 
a) # @.P eh sf 8 6 Sea OOLs Hredugel asso ops 4 . eee St ee Oe ee Pe ee {itrete's - eT ee hy Sr ee” ae ae ee Le ee a) copper rr] Py ‘ PF F ’ . 
Psd att hen halle allie Yow emrcdhla if rad rr ae ee oe re nwop OF oun a the le Oe ees OO a ee Be de MS a Ae ed ee | i qfras ’ rf ra, ran : rt ra) Fy 5 F yi A 
ar ee er ee) ot ur ew om ee er ty ey aT Oe yt Me ot at Peas ee Lah ee ey et oe oe eb ey eee ees Sry ie | ‘es eteryt oe ee ee ee ee ee ee ie ae Oy UD ae fo. rT ae . one ‘ D : 
tw Ss aoe © © Fo” ass 7 an rs Pt or ey ee ey Yor te te er et eo rYy Oe eo a ee aed i ae a be Oe Le ee nT re CGA yt Ae Tet icfos: Pt : : on : F 5 ; F A ° A Peet 
aE a OT lt See) fae we oe es Oe SMe) | A Gg Cian we a ee Lae a was n D ’ Cart ‘ - a r - ‘ Pa ' r} 
eae cab ALLA lords PCL URS oats Ue BAS UT Tal eT Av sy Oe ee ae ar ir ee Orn) Smee F Sen at o ' 
sf : CL ted ole aes | Y er | Ay en me pores c 4s ne in) tee ar) A or r fi A a a) TY 
ene rae) Py ar) om # P Pare eo ee es eee oe ir er Sor ake r ' ee) A F Sorter ate Tl ° Ar ‘ 
et © Oa meh FI OU FoF Ott Pts Ce 4 pee clip oF inc} a a e'ose prs Fa at eu? rf een Oa ‘ . o ‘a ' 
pean othe K Serle we Foa.t ets eves ye ae Oe ee ee ee ad is"ba ae oe {0 é re hr) ‘ ‘ . ee a rao 1 A A F 
. A ld wor ee te 2) Let eee LAS ee Pee, re hy et ee fey | Par} F da A Py a - A A ’ Wenn > . ‘ 
ead Not Tooke ae Pe OS Pt ee ae eet Lee a oe 2d w 8 orp hold ae a a) te rt Fi Pa = Fi . F 
‘ wate ee tL ONE ahd fe | Oe a ae et ays’, ieee ee ie oe Cet Ee tome ose 16 . 1 ey ran | Py ‘ Par = A F 
veer? ad Cel A Load Cc ie ae ere a PY ee ae) ee es va ae ee 7 Peer ae iy ny te ae ae Ok ee Oe a ae ere aon os se A ‘ ‘ ‘ae P F Ar * 
et te yang) eae i n tm Or ae Le oe oe ae ie TT ee | rin il) ee ee tee eee ee oe be ey oe Oe ee ee ee a or os * et ‘oe » 6 ’ ah . ° ’ 
flan ae cal bie oll ad de dee LTE I ice eee a ee ee ee ee ee Cr a 2 a tt en) ee | Rew ‘owe he i is Paci A aria 5 e p Py , 
ee ee ee bd , oa" 7 a ae eee) Ca a ri ree Pe re Ty ae ae ee) ee ee ee ee " F Fs “er aes A Fy ‘ * ‘ es Pa es acs 8 P ar] P 
add ee Pd oe ed et eo Pt ee ee ee te ee te et er er ne eee ver a ed go 044 Poy eae #8 gin? A eee ra a] ra a A Wie Pa 1 F rE 
aN Td kilt Le da md CA Ree ee ee ee ee ee ee ee Ls dR ool te Sad eh rea Tee Oe ee ie er ee ae Pee halt oot Pr) ny 4g ; ‘ PS oa , 5 a] 
Pie ly | P SE Se PRE ST Fant as wo G mae 8d Pe re Ye ee ae ey eT a ata) P rear - - Fi 5 Paar P A P - ; Fi : 
Lee en 7 Be ey et We ee a ere | aahesd ete oe , a Le ae a a rf ar te rere a oe se 5 Fy A F na 5 se 
ts ci) a OL TT ee eer gy) Bod te de cee¥ YY of Wy a ees ATP’ Ie har ee a nT os ae é ,? pi Prem T) Fy rie ee iar - Aa Fi 5 p Py ° 
ee Ca i a ee ee a ey Pe Pe orem lt rv) PAPC CRE TTP er ee ee a " eet ts oI P ' © ‘ p F ] . ‘ 
aa eae . aon Lee oY et ee et ee oe ee Se realy PPT RT Te tk ee rh é ri an ode tse es ' P] ort ] ] F) . 9 , ¢ F 
Lid se Ma * eet eA LY oe ee eo ee Se TL ee ee eee oe ee oe | 6 wu 1 a o.foew rd » Sep erg ss 4 Peer oe : fe Per - 5 
ode ae ye Ce tay | eee ee ea ee ieee 2 , Pate a re Or ay it Or » pels a Ce) i toe 
‘ Te ee ee et he a ee a ee ae Te eo Me oe LY >3 + a Pp ar Fi Fn A 1 
beh tt iO of Bye 0 orFeg ofeletn gp Ue heise st woe Er F ae ‘ ar) Far F Sor Py F, Py A ' 
Ue ee ee | wi wey ee ee i eee a an Pr a) PTS | F Pe ar i a OFor Pi ‘ A rare - = 7 a Ps 
Be PoE AT a a ed ee ee Se ee) Oe ee? 2 oe oe ee Oe eo Sh 4 er 5 P} ren ‘a = * = F3 ‘ *, 
Peer Me Leeda GE ct ! Oe Te a Oe ee Fra Rew sn ty Pe a | ry eo ae ’ . ORS F ar , r 4 
wt Cs Ce ee ee a ee bee Te air T A t reeves s a , AL ‘ Phar re A # A A ry 
Ce aT] em plas Ce a) p ‘ ‘ ' P Ce ir ar) 5 ° Ce ® . 
Ct re ete ee ee ee ee? a) O a OD . F ‘ ar S n A 
A) ee fob or’e na omee” rary) ae Le eee oe ee 2 oe et oe ee | ar] . A * ry r) see e P Fy 7 ry . 
PP ae ar ee Pe ee et tee ee 7 aig Tay met ey be et te ee fy F j Pa) 4 = Fs Py n oP Fy A a ‘ A 
Co ee ee | rl he ae a ae | Ly 
dos @trow od Pat LV Ee ee Sear ee Oe ay et 
roy ee eet ae ee Le ow Pee eT i Piet hy oe ee ee oe 
CIAL Le le deel Lt 









seme OLE OMe 8 dg ow tT male wae 


Fs v7 i i er ee et eee 
i a ee Ua SM ee 


















Cd ee ee 
eT La ete bee ket ee 
Pr et ee we hee 
Pree 2 Le) Fu 
i ae * Pe ae | 








© aor E 
D TP) 
tae te ae oe ee he 







roe a 
es a Pi 








a 






way 


Seren. 
Swe 









eee Fie 
bat * 
A 











e a oe tar) 

z the, 

Ce he Le oe 
.. its er 






ey vat 
Cer a) ara ts P 
ere jet -ee 


























Vv tprhe tee Ay enw a 
pees bet he Me | 68 Nip © 
hind th WT LD Peery Or TPE A DE " | i ed 









apa apy Ay Teed yt lag SL aie s Nie me vile 
Uri ee Te FL tad telah tata 
iP et) ee Te ee oe v} 
< is BL OS 4 
























. ee ed We ee td 0 he 2) ie a. re oe) tt wo gee 































































































































































Sette be Rihedealit, eta iiiclehei ded he et Sk Ae et he Se ee is arth ae ot cee ef , 4 TAN ec echt 2 . rar Sa Mir hae e A rar a p 
Ri dh ae | alii id tel bib de bes te ee TI ey Oe eh Pee Ce Ws WY he | CA ar ear | s » © Doerr) Lat er hi! rae ° . 
pepe! Read OG II Jn Sine et ied abe haada Fale LV In SN ade Ww “ert sy ya 6 : tfr,? Pi Ph a et ane A . er 5 et 
abs HO pln) io bad Medel rt oA bd aD Radind 2 Mado ¥ ea CP bbe! ‘ wees ra 4 ‘ai betes ee 4 ‘ r n f A oD ‘ ' ‘ e 
i Pi antibod Minti Peed FL ind. Wabi Bet a0 tA a dh Se be Se pas tt it at rrp et aU bee : re ere ORDEL aps Oe ae : ae ae ; ; 
ie é ve fy peddle * & Pb 1) a a] ww ea ee gt race « feats Wessel bie 6 Pa | e A Winaree s F F Sf A se nee 1 : . ” 
o . Pein a as spat Pet de abby Pind dhe he Fin 254 oily se, uy 2s ware le 0 4 a eer a ee ie Se, an) , ea Pa aan oN os 5 
O Reet ae ee LP ry eee a eet ee Sr he et Sek ee ee ran es i Al eed sree he yn ome p ar Joe i 7 ; 
ioe be te de es Be ap ere argiy are Ree yn ae) er ab. ) Per te | N - ‘ ce 5 - ae eT a ar ys an ry ; - 
til | idl Icha ela bt fod ri EEE ete weeta tie ta Sah teal lees Ti pete £ ee 4 A Snir t a 3 5 . 
a nf he ee | ae ee se oe 0 2 ee ’ ¢ ‘yw . Py B t 4 Hear f 7 ie 
Ade Ti ene hah Th tie etd bot Be a " % spon ue 7 8 : , oa : A ; Pat - 
pam Jie ave ih at Peake Oe | ela ad ple s TP P a <P 8 F a Ae . ° 
‘ {e/a Re okay pa Eb a LS SF a Te 2a a i Cr | ant “t rf an ary . r ‘ ri 5 te ry : . 
¥ Fe ratepany Here ty oy ve hs a 2 Ee ek ee 7 a Sei A ; x Ae 2 ae é Pe ; P A ane P f a : 
a he ede oe es Ee, Ld ed ee Tah Se Re ea @jVP4A ee gre . 8 ‘ "6 ‘\ . Ta nae! ee eer ' Ls H a 
« 4 a cect dy hd ee ee ed ; ey 0 4 Sere) Ya) es ey et ee eo ate arr er A P , ae ; . he ¢ H 
a Nadie bets binds LPs ELD oh Pe PLS Py Sar haere ATC oe a * Poe et ee arr te ie 5 0 “ete ae Ui era rian ae hihw ie 7 F i . 
sillerttpnd ASE alata da POO OUT LC i Lar a etre are Sybil NUT Da OR PF ie She : ° + a er) . r a a | ea D 
wa geht hic le tea ee 4 badlb-/ tlie ib} dee hy) CA, & BY seh Sid bP? i a De a OR Pe a 2 , ea ae . © ° c . Ld 
Gai Ny hi bh tie he She Th bp whined ated dk tar hh dE aE lake Me On Ah Ee, Reds ae ee ee ee) ‘ io ry, rn . A 5 oD U 
beg ly fl rare orerdp) AS ry Sed OU ed he Vase thin Sy A "erw ts aa et a : = Ae eee es P “ P P ‘ 
sod ih dielite der ae Sd | Oe ee le ee Oe Po tie? et 7) a) eee ee ee ee ae i A F - ; P . F . ls 
bie Scaled tah kd eth Wks 9 Sel Mh Mid ithe! ite td fe Be i hd Prk Aer ke ate PL he et oe Pe eT a es a) ? 3 : sf 4 ts ° ' 
A ielil Neti tal det ded edad Ta ee CT a Sie oe ee te MO RT er tT Pro) kere en oe i F a : F é ~ - 
wth Me Ts dae th ta he ET.’ eh oe er Perak eS tay 4.6 y 6 p a “ a F e 2 ee : 
hada ttn died de bi tie te ha hh Beet ee ts he La Midiad eS Sheath Titled a Mid le Mae olf Nadi Dia ttl Ad beatae | Soa Tie a ie) Fi i ‘ . ' 
. hhh ld ib Ih hel dite vin ad ee Tet ah a pap ll (Tet baa Al iia os Th aD i ae o¢ Le : ¢ i ee c 
pil ini sopaaaap det Se. 0 heey edt In el, So Moan Lon feeds Bay? ity ist eats a’ o ae) , 
<0 EES & OMETEITEPIETS ‘6 2 1b te ty he a Srey arp ena ans 0 fi o . 8 f aa 
Pee et hada ait O76 0d 08 STW mw Eg UA erty . 3 i | A ‘ ° om at 6 f . 
ape Cedi do dagen a a s + o68 rie 3 USL ed i ; 
le asco vt Bes . ne tat ey Yee ei ps o A 5 ‘ f or) 
Rime Mate de te —s r , ee 2 d . : 4 A ae ; mie 
7 Bre es tal i ee Be Tr os | OL tAh pr - .. + ae ni pn - : set 
We te othe es ee ee eve ree. ha Ve ee ’ A bot n r - a 5 
hee ee lie te Pk tt be Poe Toon eee ae ar) ha a Pa ‘ 5 D o 
ee ee ee ee ie a i eee Oe ee eA i Sa “ 2 re ; ; 
(ere ry hed pital tte wtwre mpi taw as % ew ue oe ea ake nae rea a »é oe ee Pe se ‘ + <" Pf , a ar A = aes F 
| ae ee bd deel ited WA Wd te) bad bal Mate B Dek as RITA VoOe CT A ok ie ete Th) : lige ; , aa “ ar 
Sir Prd Phen tied bbs oli Atal | w ¥ pee re > Cts ae Pe : P , o . a = . , ; 
> petal dtyese| hn. fg TOF 4 AL AL he et ee a oo ee ek ora) ee | ry 5 8 oe Ss , f v t ¥ A ny 
Aiphiith th Mrdito hata he te eth id Bell od on She ell a ota et A eR ree ae ee a) eer S ae ’ 
Lae iden AL he ee Oe te Oe Peg ee nC at ee ese re 5 PEC cn ic P un ar ae Pax ; 7 - 5 a an : 
a-Ptphhd a ee Re oe ee re earner ye uae eee Bose Da ‘- ae EE : 
bch ls tL eh ee rere beans ee mY as 4 a | | Pn 
tie eR tet Me Dy PL TV Coon © 6 Cae ’ Y Ore | a Y a 5 hl Pans o 8 
: ip lane Sth th Sadi ae Rid pt Serres Ae 64 rs ah ta Al, wheene euhsd sa se iw me : ‘ Pi ar ee ia o > 8 be ' 
phevnesdt bleh i Rin hd Atel As fin Sid Bel da he Sa) a pale a ca ea rae ™ twee 8 Fem ye ey eke He ss ec, er a] i at ie rate a Chee oa . 
"hidaed A hte Sate, tot CR BRT TL Leer i) koe Pete ee bs ee ARI I ME pls at ye ai he ie Sree eae yer eerie Gee wa ‘ . 1 ' 5, 
heel peracid By Toot Ah ok Koi UAL Ah ded oth dnt A Ok he bf a a aT Se | a) yan 1 46 Py r : 
1A, Apt > detenilly it Bete Malet TOI he an Si Ate BL AP a ee hd te ed i Bea 1 Ae Meat ary hs J AG o - ; P 
7 LBs Landihaiy i kd Zi 4 5 ' res oT aa ee Pee ery, 2 one CR 5QTD . x if i DL P ‘ a 2 c 7 





